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1.  introduction 


Difficulties  on  several  rocket  nights  in  tin;  recent  past  have  been  linked  to  the 
build-up  of  the  internal  pressure  of  the  payload.  This  increase  produced  an  en¬ 
vironment  that  was  adverse  to  the  operation  of  components  of  the  payload,  resulting 
in  tht.r  malfunction  or  failure.  To  determine  the  affects  of  these  unfavorable  con¬ 
ditions.  it  was  necessary  to  calculate  the  internal  pressure  as  a  function  of  lime. 

The  analysis  contained  herein  was  performed  to  accomplish  this  task.  It  con¬ 
sisted  of  the  mathematical  development  of  the  differential  equations  lhat  represent 
the  modelled  pressure  functions  and  the  subsequent  utilization  of  a  digital  computer 
to  determine  their  solution.  Tne  computer  programs  and  their  results  have  been 
verified  by  comparison  to  empirical  data.  These  supplementary  programs  can  be 
used  in  after-the-fact  calculations  of  the  pressure  the  payload  experienced  for  a 
particular  rocket  flight.  A  more  useful  approach  would  be  to  use  them  analytically 
to  predict  the  internal  pressure  of  a  payload  under  design.  In  this  manner,  possible 
pressure  problems  can  be  found  and  corrected  in  advance  of  the  actual  flight. 


(Received  for  publication  3  September  1980) 


2.  SINGLE  VOI  U.Wii  MATHEMATICAL  ANALYSIS 


The  problems  from  pressure  build-up  arise  from  the  fact  that  the  internal  gas 
cannot  vent  fast  enough  to  lower  tne  payload  pressure.  A  pressure  differential  is 
formed  which  results  in  loads  to  sensitive  items  such  as  doors,  compartments  and 
components,  causing  damage,  malfunction,  or  inoperation  of  these  devices. 

On  the  launch  pnd,  the  payload  is  usually  pressurized  to  insure  that  the  "clean" 
area  is  at  a  slightly  higher  pressure  than  the  surrounding  environment.  In  this 
manner,  dust  particles  can  be  kept  from  contaminating  the  important  payload  areas. 
As  the  rocket  ascends  through  the  atmosphere  the  external  pressure  drops  faster 
than  the  payload  pressure  can  follow,  creating  a  differential  between  the  internal  and 
external  pressures. 

The  payload  is  modelled  as  a  simple  box  of  volume  equal  to  the  volume  of  gas  to 
be  vented.  To  this  box  is  attached  the  venting  apparatus  of  the  payload,  as  shown 
in  Figure  1.  The  configuration  consists  of  any  combination  of  valves,  filters  and 
orifices  used  in  relieving  the  internal  pressure  o'  the  payload.  The  only  way  for  the 
internal  pressure  of  the  payload  to  be  released  is  through  the  venting  apparatus  and 
any  leaks.  With  the  proper  design  configuration,  the  internal  pressure  can  be  made 
to  approximately  track  the  external  pressure. 

The  valves  are  closed  until  the  pressure  differential  across  them  is  equal  to  or 
greater  than  their  cracking  pressure.  The  filters  are  open  and  operating  at  all 
times.  Provision  is  included  for  any  leaks  due  to  doors  or  seals,  which  are  modelled 
as  orifices  with  an  effective  exit  area.  More  will  be  presented  on  the  operating 
characteristics  of  these  devices  later  in  this  report. 


Figure  1.  Payload  Model 
for  Venting  Analysis 
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The  mathematical  analysis  is  begun  by  writing  the  perfect  gas  law  for  the  fluid 


in  the  bos: 


p.  =  p.RT. 

i  i 


m.rfT  /V  . 
1  i 


Differentiating  with  respect  to  time,  we  obtain 

dp.  ,  m  RT  p  dm.  dT 

<IT  I  - V—  =  V  Ti3T-  +mi  dt“  l2) 

L  J  L 

For  simplicity,  the  internal  temperature  will  be  assumed  be  constant  throughout 
the  flight.  The  average  ascent  through  the  atmospher  lasts  80  sec;  this  does  not 
allow  enough  heat  to  be  transferred  to  or  from  the  gas  to  significantly  change  its 
temperature.  Then  Kq.  (2)  becomes 


RT.  dm. 
s  i 


where  dp./dt  is  negative  for  p.  >  p(  -  To  calculate  the  internal  pressure  from  Eq. 
(3),  we  need  to  determine  the  rate  of  change  of  the  mass  in  the  box.  The  fluid 
mass  is  decreased  by  tiie  flew  through  the  venting  apparatus,  which  is  governed  by 
the  flow  characteristics  of  the  valves,  filters  and  orifices.  These  characteristics 
have  been  measured  experimentally  from  actual  hardware  and  arc  presented  in 
detail  in  Appendix  A. 

The  mass  flow  rate  is  a  function  of  the  density  of  the  gas  and  the  pressure 
differential  across  the  oeviee,  which  is  defined  as: 

AP  -  -’’int  -  pext  '  (4) 

The  external  pressure  is  taken  as  the  atmospheric  pressure  (at  altitude)  that  the 
payload  experiences,  which  can  be  determined  from  the  rocket's  trajectory. 
Aerodynamic  affects  producing  a  pressure  coefficient  and  a  subsequent  change  in 
the  "external"  pressure  are  ignored  for  the  following  reasons: 

(1)  In  all  cases  analyzed  so  far,  the  valves  and  filters  have  been 
mounted  on  the  cylindrical  sections  of  the  payloads,  for  which 
Cp  is  negligible  or  zero: 

(2)  The  effects  of  angle  of  attack  and  boundary  layc"s  are  considered 
to  be  negligible. 


The  total  mass  How  rate  is  the  sum  of  the  contributions  of  the  individual  venting 
components: 

=  n  (b  of  valves)  +  m,  (#  of  filters)  m  .  (S) 

IV  f  o 

From  the  considerations  of  continuity,  we  find  that  the  rate  of  change  of  the  internal 
gas  mass  is  equal  tc  the  rate  at  which  the  gas  leaves  the  volume: 


mi  =  ~mT  * 


Using  this  substitution,  the  differential  equation  which  models  the  payload  internal 
pressure  becomes 

do.  RT. 

1  i  _  1  - 


With  the  external  pressure  as  a  function  of  time  and  the  equation  for  the  mass  flow 
rate,  the  above  differential  equation  can  be  integrated  to  give  the  internal  pressure 
as  a  function  of  time. 


:t.  COMPUTER  PROGRAM  PRESSt.FOR 

Equation  (7)  can  be  integrated  numerically  to  obtain  the  internal  pressure  at 
any  through  the  use  of  a  digital  computer.  The  general  ordinary  differential 
equation  of  the  form  dy/dx  =  f(x.v)  with  initial  condition  y^  -  U*Q).  is  solved  using 
a  fourth-order  Runge-Kutta  integration  process.  This  is  a  single-step  method  in 
which  the  value  of  v  at  s=x  is  used  to  compute  y  ,  a  v(x  ,).  The  relevant 
formulas  for  integration  are: 


=  v  *■  —  (T  -  Tf  t-  2T  *  T  ) 
-n  6  1  2  3  4 


hf(x  ,  v  > 
r  '  n 

hflx  +h/2.  v  +T./2) 
n  -  n  1 

hf(x  -rh/2.  y4T.}/2) 
n  n  2 

ilf(x  +h,  v 

n  -  n  4 

step  size  , 


A  step  size  of  0.  0»  hi, c  tvus  selected  for  this  program.  The  schematic  of  Figure  2 
gives  a  visual  representation  of  the  integration  process.  The  y  values  defined  at 
intermediate  steps  are  used  to  calculate  the  Suture  y  values,  without  using  any 
previous  results.  The  method  is  a  quick  sac  efficient  means  of  Integrating  . 
differential  equation. 

This  integration  technique  raid  the  oiffenential  equation  Save  been  programmed 
on  *he  PDP-11/34  computer  using  the  FORTRAN  IV  language.  A  listing  xhd  flow 
cii'i rt  of  program  PRKSS4.FOR  are  presented  in  Appendix  B. 

X  DOMAIN 
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Figure  2.  integration  Schematic- 
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4.  MULTIPLE  VOLUMES  MATHEMATICAL  ANALYSIS 


[  i 

II  1. 


A  payload  can  be  constructed  such  that  interconnected  compartments  are  in¬ 
volved  in  the  venting  process.  The  model  Of  this  case  consists  of  boxes  connected 
by  various  venting  components,  as  shown  i r.  Figure  3  for  two  volumes.  The  main 
volume,  box  No.  1,  is  set  up  in  a  manner  similar  to  the  single  volume  model: 
with  valves,  filters,  and  orifices  exposed  to  *he  external  environment.  The  major 
difference  is  that  it  now  has  fluid  input  from  *hc  secondary  volume.  The  secoodaiy 
volume,  box  No.  2,  is  sat  up  for  venting  gas  to  both  the  main  volume  {through 
valves)  and  the  external  environment  (through  filters  and  orifices). 

For  volume  No.  1.  the  mass  flow  rate  is 


m.  =  — (m_  -  m_  ) 

*  1  1  *  O 


where  mT  is  tl  '1  flow  rate  output  of  the  main  venting  apparatus  and  m_,  is  the 
1 1  l2 
rate  input  from  the  second  volume.  For  volume  No,  2,  the  mass  flow  rate  is 


if  p 


Proceeding  in  a  namwr  analogous  ic  the  single  volume  dtalysis.  we  obtain 
differentia!  equations:  o.h-  for  each  volume- 


Thi  alysis  can  be  extended  to  include  an  infinite  number  of  volumes,  obtaining 
n  differential  equations  for  n  volumes  interconnected  together  and  venting  to  the 
atmosphere. 


5.  COMPUTER  PROGRAM  PRESSM.FOR 

A  two-volume  venting  calculation  has  also  been  programmed  on  the  PDP-11/J4 
computer.  The  numerical  integration  process  is  similar  to  that  used  for  the  single¬ 
volume  problem;  here  it  is  set  up  for  a  system  of  ordinary  differential  equations. 
The  method  is  contained  in  the  FORTRAN  IV  subroutine  RKGS.  FOR  in  the  IBM 
Scientific  Subroutines  Manual  (Reference  4).  A  lit  'g  and  flowchart  of  program 
PRliSSM.  FOR  arc  presented  in  Appendix  C. 


6.  SAMPLE  PROGRAM  RESULTS 

The  two  computer  programs  referenced  in  this  report  have  been  used  to  ana¬ 
lyze  payload  designs.  Presented  here  are  the  results  for  the  SPICE,  IRBS,  and 
ZIP  payloads.  Also  presented  are  the  results  of  the  test  case  used  to  verify  pro¬ 
gram  operation  and  validity. 

6.1  SPICE  Payload 

The  SPICE  payload  was  flown  on  27  Jan  1979  and  experienced  failure  of  the 
door  unlatching  mechanism.  The  cause  of  this  failure  was  determined  to  be  a  build- 
uo  of  the  internal  pressure  which  resulted  in  increased  loading  on  the  door.  The 
SPICE  payload  was  analyzed  using  an  earlier  version  of  the  single  volume  computer 
program.  Figure  4  is  a  plot  of  the  computer  results. 

SPICE  payload  configuration: 

Volume:  19.55  cu  ft. 

Venting  apparatus:  2  P7-637  0.  50  psi  relief  valves. 
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6.2  IRBS  Payload 

The  IRBS  payload  was  analyzed  during  the  testing  phase.  It  consists  of  a  small 
chamber  venting  into  the  larger  main  volume  of  the  payload.  The  intent  here  was 
to  determine  the  maximum  internal  pressure  that  the  payload  would  experience  in 
order  to  generate  proper  testing  levels.  Figure  5  is  a  plot  of  the  program  results 
for  Volume  No.  1;  Figure  6  shows  those  for  volume  No.  2.  Notice  that  the  valves 
on  the  secondary  volume  do  not  operate  at  all;  this  is  shown  bv  the  secondary  volume 
internal  pressure  being  lower  than  the  main  volume  pressure  in  Figure  6.  The  leak 
takes  care  of  any  pressure  build-up  in  volume  No.  2. 
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PRE55URE  -  <P5I> 


IRBS  configuration: 

Volume  No.  1: 
Venting  apparatus: 

Volume  No.  2: 
Venting  apparatus: 


46.  80  cu  ft, 

3  P7-637  0.  50  psi  relief  valves  venting 
to  the  atmosphere, 

0.  18  cu  ft. 

2  P-249  0.  10  psi  relief  valves  venting  to 
the  main  volume;  1  leak  (orifice) 
with  effective  area  of  0.  000042  sq  ft 
venting  to  the  atmosphere. 
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1RB5  PRYLDRD 
SECOND Y  VDLUME 


Figure  6.  Program  Results  for  IRHS  Payload-Secondary  Volume 
6.3  /IP  Pa\loail 

The  ZIP  payload  was  analyzed  during  the  design  phase.  It  is  a  good  example 
of  how  the  programs  can  be  used  to  pin-point  problems  in  advance.  The  original 
analysis  of  ZIP  showed  a  maximum  delta  p  of  over  1. 69  psi;  it  also  showed 
0.  187  psi  at  the  critical  time  of  door  unlatching  and  opening  (see  Figure  7).  The 
venting  configuration  was  revised  to  include  the  larger  relief  valves.  Subsequent 
leunaHsis  showed  a  much  improved  situation;  the  maximum  delta  p  was  decreased 
to  .just  over  0.  50  psi.  Figure  8  presents  the  revised  pressure  prediction  for  the 
ZIP  payload. 


ZIP  configuration: 

Volume: 

6.  80  cu  ft. 

Original  venting: 

6  P-249  0.  10  psi  relief  valves, 

2  RA-2500  filters  0.  11045  sq  in.  in  area, 

Revised  venting: 

10  P-249  0.  10  psi  relief  valves, 

4  P7-637  0.50  psi  relief  valves, 

2  RA-2500  fi’ters  0.  11045  sq  in.  in  area. 

ZIP  PRYLDRD 
□RIEINRL  CONF 

NS\\  IWTERNHL 

BELTS  P  — 

X 

TIME  -  <5E 


KM 


Figure  8.  Program  Results  for  ZIP  Payload— Revised  Venting  C.'onfiguration 


6.4  Program  Test  Configuration 

In  order  to  validate  the  pressure  programs,  a  test  case  was  developed  and 
evaluated  during  the  ZIP  payload  analysis.  A  standard  volume  with  only  one  relief 
valve  was  evacuated  using  a  vacuum  pump,  such  that  the  external  and  internal 
pressures  were  known  to  an  accuracy  of  ±  0.  1  psi.  The  tost  case  results  -\i e re  the: 
compared  to  those  predicted  by  the  computer  program  for  the  same  external  pres¬ 
sure  variation.  This  comparison  lead  to  additional  refinements  in  the  programs, 
with  subsequent  improvement  in  their  prediction  capability.  A  comparison  of  the 
results  follows  in  Figure  9.  It  is  evident  that  the  program  has  sufficient  accuracy 
for  design  work  while  remaining  slightly  conservative. 

Test  configuration: 

Volume:  1.  00  cu  ft. 

Venting  apparatus:  1  P-249  0.  10  psi  relief  valve. 
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Appendix  A 

Operating  Characteristics  of  Venting  Components 

The  venting  apparatus  is  analyzed  as  an  opening  of  a  certain  area  through  which 
the  fluid  flows.  The  mass  flow  rate  takes  the  form  of 

m  pQ  5  pvA  .  (Al) 

The  density  of  the  fluid  is  a  function  of  the  pressure  and  the  temperature  and  will 
be  determined  by  inlet  and  outlet  conditions.  The  volume  flow  rate  depends  upon 
the  velocity  of  the  fluid  and  the  area  of  the  opening.  Both  of  these  are  functions  of 
the  fluid  pressure.  The  velocity  is  related  to  the  pressure  ratio  across  the  open¬ 
ing.  In  the  case  of  the  relief  valves,  the  opening  area  is  variable  and  is  dependent 
upon  the  pressure  differential.  The  cracking  pressure  on  a  valve  is  controlled  by 
the  properties  of  the  helical  spring  which  is  part  of  the  valve  mechanism  (see 
Figure  Al).  As  the  delta  p  is  increased,  the  valve  opens  and  the  area  of  the  open¬ 
ing  is  dependent  upon  how  the  spring  is  compressed.  At  some  value  of  the  pressure 
differential,  the  valve  will  "bottom  out;"  that  is,  the  spring  will  reach  maximum 
compression  and  the  exit  area  will  be  at  its  greatest  value. 

We  must  also  take  into  account  compressibility  affects  and  the  phenomena  of 
choking  when  analyzing  the  venting  apparatus.  The  flow  through  an  orifice  (or  any 
opening)  can  increase  its  velocity  only  until  the  Mach  number  reaches  the  value  of  1, 
At  this  point  the  velocity  in  the  throat  (smallest  area  of  the  orifice)  becomes  sonic 
and  the  volume  flow  rate  reaches  its  maximum  value.  Any  attempt  to  further 


21 


increase  the  velocity  through  an  increase  in  della  p  will  not  be  successful  in 
changing  the  volume  flow  rate. 


Pi  Pi 


SPRING 


Pe 

f*"  VALVE  HEAD 
O-RING 

MOVE  SEAT 
VtttyE  BODY 


Figure  A 1.  Throat 
Concept  and 
Valve  Mechanism 


This  is  so  because  an  increase  in  velocity  would  cause  the  Mach  number  at 
the  throat  to  be  greater  than  1,  which  is  an  impossibility  from  fluid  dynamics. 
When  the  velocity  in  the  throat  equals  the  speed  of  sound,  the  flow  is  said  to  be 
choked.  The  volume  flow  rate  will  remain  constant,  no  matter  how  large  the 
pressure  differential  becomes. 

Choking  is  a  function  of  the  pressure  ratio  across  the  opening;  here  pe/p..  As 
the  velocity  increases  and  M  approaches  1,  the  pressure  ratio  decreases.  At  a 
certain  value  of  pe'p.  the  flow  will  reach  M  1;  this  is  railed  the  critical  pressure 
ratio  and  signals  the  onset  of  choking.  A  numerical  value  for  the  critical  pressure 
ratio  can  be  determined  from  fluid  dynamics. 

Writing  Bernoulli's  equation  for  compressible  flow,  we  have 

o  9 
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The  velocity  well  away  from  the  opening  inside  the  volume  is  negligible;  thus 


v.  =  0  and 
1 


_  dp.  v"  „  dp 

I  _li  i  i  f  _UL 
J  P;  P„  * 


Rearranging  the  above  equation,  it  becomes 


v"  ,  do 

_ e  _  r  dp  ,  f  e 


e  -  /HE;  _  f  o 
9  J  O-  J  P  ' 


}J  f 


From  the  relations  of  isentropic  flow  (frictionless,  adiabatic  flow  of  a  perfect  gas) 
we  have  the  following: 


P  =  Cp‘  p=  (p  lC)lly  C 


Uy  c  1  ly  -  £ 

P 


Using  this  in  Eq.  (A4),  it  becomes 


(1  ! 

i  :!  i 


s 


)j  1 


-2s-  =  J(C/Pi)1/ydP.  -  f(c/pe)1/y  dpe 

1  .  CI/>'  (/dp,/^  -  /-P./P^) 

y”  P,l/»  ( P“'1,V>  p“-l/y>’S 

2  e=  Pj  y  '(T-'i'/r )'  (1-1 /y )'  y 


%Ur*lr- 1)  (Pri,/r  -  Pri)/r 


!_S.  =  y  Uy -l)  1L  f(l-(pe/p.){>"1,/r 


Therefore,  the  exit  velocity  of  the  opening  is  given  by 


2  v  D;  / 


v0  \  TT  p, 


At  Mach  =  1,  v  will  equal  the  speed  of  sound  and  p  Ip.  =  (p  /p.)  then: 

c  Cl  el  crn 


v  =  a  = 
e 


(l 

1  Pe  t  r-i  Pi  ^ 


i-(p  lp.)ir'l)lr 
*e  i cnt 


(A  12) 


Rearranging  Eq.  (A  12)  gives: 

i  -  <p /p.)(r:I)/>'  =  V-  ^  ^ 

C  i  crit  ^  n  n. 


2  pc  Pi 


(A  13) 


but 


£  =  <Pi/-°c)1/r 


(A  14) 


from  the  isentropic  relations  used  earlier.  Substituting  Eq.  (A  14)  into  Eq.  (A13) 


1  _  (p  /p.)<>'  -  »>/> 

*e  *  i  cnt 


r-i 

2 


V^crit 

(p/p.)1  ^ 
c  l  cnt 


<Pe/pi)crit 


D/r 


(A  15) 


(l-f(y-l)/2)(pe/p.)^r:t1)/>'  =  1  (A16) 

and  solving  for  the  pressure  ratio,  we  have 

<p  Ip.)  ..  =  <2/(y  +  l))''/(v_1)  .  (A  17) 

'  e  I  cnt 

Equation  (A.-)  is  the  expression  for  the  critical  pressure  ratio.  For  a  y  -  1.40, 
the  critical  pressure  ratio  is  0.  5283.  When  the  pressure  ratio  is  less  than 
(p  Ip.)  the  flow  is  choked  and  the  volume  flow  rate  is  maximum.  Eqs.  (All) 
and  (A  17)  have  been  used  in  the  computer  programs  referenced  in  this  report. 

The  volume  flow  rate  for  the  venting  apparatus  used  on  our  payloads  has  been 
measured  experimentally.  The  valves,  filters,  and  orifices  are  modelled  by 
developing  mathematical  formulas  for  their  flow  characteristics.  Originally,  this 
was  accomplished  from  manufacturer's  data.  However,  upon  close  examination 
that  information  was  found  to  be  quite  dated  and  it  was  deemed  necessary  to  expert 
mentally  test  the  apparatus.  The  numerical  information  on  the  valves  and  filters 
presented  in  Table  A1  is  based  upon  this  empirical  data.  Leaks  are  modelled  as 
orifices  and  arc  governed  by  the  following  equation. 
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m  =  A  ^64.  S/)Ap 


F rom  which  we  see  that  Q  -«  Ap,  or  that  InQ  «=  0, 5  2nAp  for  an  orifice. 


Table  Al.  Valve  and  Filter  Mathematical  Models 
CIRCLE  SEAL  Pressure  Relief  Valves 

James.  Pond  &  Clark,  Inc.  s  Pasadena,  CA _ _ _ 

P-24S-U. 10 


Labeled  cracking  pressure- 

0.  id 

psi 

Measured  cracking  pressure 

0.  0387 

psi 

Knee  pressure 

0. 10 

psi 

Volume  flow  rate  curve: 

Ap<  0.  10  (nQ  =  10.8789 


Apk  0.  10 
-0.  50 


0.9707  ~  0. 4930  InAp 


Labeled  crackii:g  pressure  0.50  psi 

Measured  cracking  pressure  0.3250  psi 

Knee  pressure  0. 59  psi 

Volume  flow  rate  curve: 

Ap<  0.  59  fnQ  =  12.7900  +  17.  3978  fnAp 
Apk  0.59  -  3.8547  -*•  0.4783  fnAp 

Valve  vciume  flow  rate  modelled  by; 

Q  =  exp(A  *  R  (r*Ap) 

TmILLIFORE  Filters 

Milliporc  Corn. ,  Bedford,  MA _ ____ 

CIV- 19  Cartridge  Filter 
A  =  —  5. 6347 
D  =  114.9396 
C  =  -60.  04  16 
D  =  12. 9680 

Filter  volume  flow  rate  mo  felled  by: 

Q=A  +  BiAp)  +  C(Ap>2  -  OlAp)3 

Filter  length  Q  multiplier 

31  inch  Q  X  1.4! 

22  ir.ch  Q  X  1.  00 

12  inch  '  QX  P.50 

RA-2500  Membrane  filter  1.2  pm  pore- size 

A  =  -0. 00702 
B*-.  2.0191 

A  =  Exit  area  in  square  inches 

Filter  volume  flow  rate  modelled  by: 

Q  =  [A+B(Ap)JA 


.  The  knee  preSSure  Usted  *  Table  Al  is  that  pressure  at  which  the  valve 
oottoms  out  and  its  flow  area  becomes  constant.  The  valve  will  then  begin  to  act 
m  a  manner  similar  to  an  orifice.  This  is  reflected  in  the  values  of  the  slopes  of 
e  flow  curves  approsunately  equalling  0.  5.  Any  discrepancies  a„  probablv  due 
to  ‘he  ejects  of  discha-e  coetficients  which  are  not  directly  taken  into  account 
here.  They  are  a  function  of  the  fluid  pressure  and  the  Reynolds  number. 

the  f3Ct0r  ^  ^  3PPUed  t0  thC'  3b°Ve  natBerical  data  to  correct  for 

tne  lact  that  the  measu rements  were  -  t  -  *  - 

ere  ^en  atmospnenc  pressure  and  are  bem® 

applied  at  altitude  (lower  than  atntosphsric  -'-essurel  e 

F  r^saurei.  r  rom  the  orifice  Eo.  CA1SI 

we  see  that  Q  is  mainly  a  function  of  the  density; 


Q  =  f(y  pAp  Ip) 


(A  19) 


Rearranging  we  have 


Q  =  f(p1/2/p)  =  f(l/pi/z)  =  f(i/p1/2) 


<A20> 


Thus*  the  volume  flow 
pressure: 


rate  is  inversely  proportional  to  the  square  «Mt  of  the 


Qatmai'Pato  Hpj/2 


IA2V, 


Then 


Q./Q 

V^aijn 


=  {patm/pi> 


1/2 


(A  22) 


«md  the  low  pressure  correction  is 


Qj  =  Q_.  (p  . 

»  atm  ratm 


/p.-> 


1/2 


(Alai 


where  Qatm  13  the  measured  volume  flow  rate  at  atmospheric  conditions.  The  i» 
pressure  correction  was  included  with  the  compressible  flow  conations  in  the 
computer  programs. 

A  comparison  of  the  computer  prediction  and  the  measured  «sute  e,  ^  teJt 
configuration  .boned  tta.  tbe  pregrem  ^  «C5p,  „  ^  ^ 

(lower  cflcreai  preuuurcul.  Eiperimeufing  „*  computer  ^  ^ 

ound  that  an  additional  correction  factor  of  the  form 


Q.  —  Q (p./p  )n 

i  atm  *VFatm' 


.’Will . . .  .i>il|l'!!l!>|l'l|i|[l||{||i|i!lililii|llill|i|l|li|l||ii'll' . . 


(where  n  is  between  0.  20  and  0.  30)  increases  its  accuracy.  This  additional  correc¬ 
tion  factor  could  possibly  incorporate  the  affects  of  a  discharge  coefficient.  The 
total  correction  factor  applied  to  the  measured  volume  flow  rate  is: 

Qcorr  =  50  “  *  (A25> 

The  following  then,  are  the  operating  conditions  that  are  used  in  the  programs 
developed  for  the  PDP- 11/34  computer. 

For  valves: 


Ap  <  pc 

Ap  2  p  p  / p.  >  (p  /p.)  .. 

*  *c  *e  *i  *e  *i  ent 


Pc/Pi^Pe^crit 


For  filters: 


Al  <  1 


lM  =  1 


Q  =  f(Ap) 

Q  =  const. mt 
(choked) 


i\l  <  1 
Q  =  f(Ap) 

M  =  1 

Q  =  constant 
(choked) 


m.  =  pQQ 
f  ^  corr 


m  0 
v 

m  =  p  QQ 
v  ‘^corr 


m  =  pQ 
v  pn 


previous 


nif  =  PQpre 


For  leaks  and  orifices: 


i\l  <  1 


M  =  1 


Q  =  f(Ap) 

Q  =  constant 
(choked) 


previous 


m  =  pQ 
o  r 

m  =  pQ 
o  '  'previous 


From  the  above,  we  can  see  that  even  though  the  venting  apparatus  becomes  choked, 
the  MASS  flow  rates  can  increase  or  decrease  because  the  density  of  the  flow  can 
increase  or  decrease.  Only  the  VOLUME  flow  rate  is  affected  by  choking. 
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Appendix  B 

Listing  and  Flowchart  of  PRESS4.FOR 


Computer  Program  PRESS4.FOR 


Language;. 

Computer; 

Memory  Requirements:- 
Port  ran  File  Size: 


PORT  RAN  IV 
DEC  PDP-11/34A 
8K  Words 
2 1  Blocks 


Input  Pile:- 
Output  Files: 


PTN30.  DAT 
PTN31.  DAT 
FTN32.  DAT 


Major  liquations  Used:-  Equations  (7),  (All),  (A18), 

(A25),  and  those  of  Table  Al. 

Integration  Technique:  Fourth-order  Runge-Kutta 

process,  Eqs.  (8)  and  (9); 
see  also  Subroutine  RK2 
of  Reference  4. 


« 
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PROGRAM  PRESS4 
I.' - PROGRAM  F'REHS4 .  F  OR 

C  --  -IHFS  PROGRAM  IS  AN  AFFFMPf  TO  PRFDICI  FHE  tN I ERNAL  PRESSURE  OF  A 
C-  — PAY! UAH  Of  A  kOCKR!  AS  ITt  ASCENDS  THROUGH  1  HE  A 1 MOSPHERC .  I  HE 

L - PAYLOAD  IS  SEI  UP  AS  A  CHAMBER  WITH  ATTACHED  VENTING  VALVES. 

C - PILFERS  AND  ORIFICES.  HOUR  LEAKS  ARE  MFJDELF-ED  AS  ORIFICES.  FHE 

C - PROORAH  NUMERICAL l.Y  IN1ECRAFFS  1  HL  11 CFFEKEN  I  X Al.  EOFJA *  l(JN  ( HA  I 

C - REPRLSfc N  F  S  I  HE  PRESSURE  HER  I VATl VE . 

C - - — AUFHUR!  C.  P,  KRESS 

Xi  I  MENS  I  UN  GAS (5)  .  I  I  It  £<20>  .PRLXI  <  JOO)  .PF  CMC  (  IOO) 

RE  A I  MASS. MULT 

DA  F  A  ORA  XF  1  .  ORA  i£2 .  ORA  TEE' ,  (FRA  I F  L  »  RA  IF  *  RA  t  EF  » RA  f  F  F  .  RA  I F-  V  F .  PA  I  fc  V2  i 
Y  DCUEFF .  INDEX. I  INDEX. I  MARK/9Y0.  . 0 . VO . 2 , 0 r 0/ 

l, - FORMA  I  H I A  FEME N 1 S 

9000  FORMA!  <'»0A4) 

91  00  FORMAT f I  ‘0./F10.0) 


9 1 v |  1  ORMA  F  > 

(F  10.0' 

9200  FORMA  Ft 

'10X. 

'HAS  PROPER F FES i '/ 

* 

UlX. 

1  I PF.  : 

' .bA4/ 

Y 

1  5X  . 

MAIN  VOLUME 

'  - 1-  F  0 . 2 ,  ' 

('ll  FI'' 

Y 

FSXf 

'INF  UAL  PRESSURE 

= 

'  -F  10.2.  ' 

T'S  t '  / 

Y 

ISX. 

TEMPERA FURE 

'  »  F  1 0 .,  1  r  ' 

DEGREES  !  '/ 

Y 

FSX. 

GAS  CONS r AN F 

' .F  10. 2« 

F  (  -1  B/LB-DEG 

Y 

//) 

9201  FORMA! </IOX? 

VAF  vf  onf:  PROPER  !  I  ES : '  / 

Y 

I5X. 

i  ypl: 

'  •  2‘>A4 ) 

9202  FORMA  1 l 

1  5X  > 

'NUMBER  OF  RELIEF  VALVES 

= 

' • FiO . 0/ 

Y 

15Xr 

'CRACKING  PRESSURE 

= 

' .FIO. 2* ' 

PST'/ 

* 

F5X  r 

'CURVE  CHANGE  POINT 

- 

'  .  F  1 0 . 2  ?  ' 

PSI'/ 

Y 

1 5X  > 

'rOEPFICIF„Nf  1 

- 

'  .F  F0.3/ 

Y 

FSX. 

' COEFF 1CXEN F  2 

'  .  F  FO  .  3/ 

Y 

1SX. 

COEFFICIENT  3 

- 

'  •  FlO . 3/ 

Y 

15X. 

'COEFF 1CIENF  4 

•= 

' .FIO. 3/) 

9203  FORMA F < /l OX. 

•valve  two  properties: '/ 

* 

ISXr 

'  1  YPE : 

' 7 20A4 ) 

9206  FORMAT </TOXr 

'FILTER  PR0PER1 TES: '/ 

Y 

lSXr 

'  I  YPE ! 

' .  20A4 ) 

9207  FORMA T< 

15Xf 

'NUMBER  OF  FILTERS 

' .FI 0.0/ 

Y 

15Xf 

'EXIF  AREA 

= 

' rFlO.S,  ' 

SO  IN'/ 

Y 

tSX. 

'COEFFICIENT  1 

- 

' . FIO . 3/ 

Y 

ISXf 

'COEFFICIENF  2 

- 

' .FIO. 3/ 

Y 

lfvX  F 

'COEFFICIENT  3 

= 

' .FIO. 3/ 

Y 

1 5X  f 

'COEFFICIENT  4 

' .FIO. 3/) 

9209  FORMA T</1 OX f 

'DOOR  LEAK  PROPERTIES: '/ 

Y 

15X  f 

'TYPE  OF  SFAI.: 

'  .  20A4 ) 

9210  FORMA  IF 

XSXf 

'EFFECTIVE  AKtA 

' .FIO. 8.  ' 

SO  FI '/) 

9212  FURMAI <//IOX 

. 'CHOKING  PROPER! TES! V 

Y 

TSX 

. ' RATIO  OF  SPECIFIC  HLA I 

s 

=  '.HO.  3/ 

Y 

F5X 

i  '  CRI F ICAL  F-RESSURE  RA  1 10 

=  '.FIO. 4/ 

Y 

15X 

.  ' SPEED  OF  SOUND 

=  '  . F 10, 1 

.  '•  FPS'///) 

92 1 5  FORMA!  U2X. 'EXTERNAL  f  N  TERNAl  ' .  1  4X . ' TNT  ERNAL ' .  4X  • ' TOTAL  MASS'/ 

*  SX.'TIME  '.2('  PRESSURE' > »3X. 'DELTA  P',4X.'fiAS  MASS ' . 

Y  SX-'FLOU  RATE' /5X. 'SECS '  .6X .  'PSI ' »2<  7X. 'PSI  '^  . VX. 'I  DM' . 

*  8X. 'LBM/SEC'/) 

9216  FORMAT (K9.1r2F10.2.F10.3.F13.SjF14.7) 

9220  F0RMATIS6X. 'VALVE  ONE'- .9X. 'VALVE  TWO '  .OX  .  '  F  XL  I  ER  ONE  '■  . 

%  9X.LEAK  0NE'/13X.'P  R  F.  S  S  U  R  E ' » 1 9X » 4  <  9X  »  '.FF.OW  RATE')/ 

Y  2X  f  'TIME  EXI  INI  RAF  10  D1FFR' .3X. 'MACH  IiENSTTY  '• 

*  4< 'VOLUME  MASS  ' )/2X • 'SECS '  *2(3X> 'PSI ' >  »3X» 'PE/PI ' . 

*  '  PS  I  NO  LDM/C1J  FT  '.4('  CF/S  LBM/SEC  ')/) 

9221  FORMAT  <F6.1fTXf2F6.2fF7.4fF7.3fF8.3.F9,5f4<F8.4fF10.6) ) 
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n  o 


0 - 1 1  HE  REFERENCES 

1 IME-O.OO 

c: - ALLOCATE  DATA  FILES 

N1N=30 

NOUI-NIN+1 

N0UT2=N0UT+1 

C - INPUT  AND  OUTPUT  OF  VITAL  INFORMATION 

C - GAS  PROPERTIES 

READ (NINr 9000)  TITLE 
READ<NINr9000)  GAS 

READ(NINr9101 )  VOLUME r VOL r PINT r  TEMP r RGAS r GAMMA 
WR I TE ( NOUT  r  9000  >  TITLE 

WRI TECNOUT f 9200)  GASrVOLUMErPINI r TEMP r RGAS 
WRI IE(NOUI 2 r 9000)  TITLE 
WRI IE(N0UT2r9220) 

0, - VALVE  ONE  PROPERTIES 

50  READ (NINr 9000)  TITLE 

READ (NINr 9101 >  VALVS1 r CRACK 1 yCHANGl rAl r A2r A3 • A4 
WRI TE ( NOU I r  9201 )  TITLE 
IF  (VALVSt.EQ.O. >  GO  TO  100 

WRIT E ( NOU Tr 9202)  VALVS1 r CRACK t  r CHANG! .A1 rA2r A3 r  A4 

C - VAt  VE  TWO  PROPERTIES 

100  READ (NINr 9000)  TITLE 

READ(NINr9101 )  VALVS2 r CRACK2 r CHANG2 r H 1 .B2rB3rfa4 
WRlTE(NOUT r9203)  TITLE 
IF  (VALVS2.E0.0. )  GO  TO  1 1 0 

WRIT! (NOUT  r 9202 >  VALVS2 r CRACK2 r CHANG2 r B1 r H? r B3r HA 

C - 1-  il  TER  ONE  PROPERTIES 

1J0  READ (NINr 9000)  TITLE 

READ(NINr9101 )  FIL  TRS  r AFILT .Cl.C2rC3fC4 

WRI TE(N0UTr9206)  TITLE 

CF  (F1LIRS.E0.0. )  GO  TO  120 

WR  1 1  E  ( NOU  T  r  9207 )  FIL  I RS  r  AF I L  T  r  f,  1  r  C2  r  C3  r  C4 

C - LEAK  ONE  PROPERTIES 

120  READ (NINr 9000)  TITLE 

READ (NINr 91 01)  ALEAK 
WRI TE(NOUTr 9209)  TITLE 
IF  ( ALEAK. EQ. 0. )  GO  10,170 
WRITE(NQUT ?9210)  ALEAK 

C - INPUT  OF  EXTERNAL  PRESSURE  HISTORY 

170  READ ( NINr  9101 )  1 TMENDr TS TEP r DIVIDE 

READ ( NIN r 9100 )  MAGNI r 

RE AD ( N I N  r  9 1 0 1 )  <  PTIM£( N  >  r PREXT (N  >  r N=1 rMAGNIT) 

C - CONVERSION  TO  PROPER  UNITS 

TEMP=  TEMP+459 . 67 
PINT=PINT+PREXT(1 )/144. 

PA  TM=PREXT < 1 ) 

CALCULATION  OF  INITIAL  CONDI  TONS 
INTERNAL  GAS  DENSITY 
MULT=RGAS#TEMP/ VOLUME 
MASS=PINT*144./MULT 
DENST  Y=MASS/VOLUME 

C - PRESSURE  DIFFERENTIAL  AND  PRESSURE  RATIO 

PEXT  =PREXT ( 1 ) /144 . 

DELT  A=PINT-PEXT 
PRA  TIO=PEXT/PINT 

c - CRITICAL  PRESSURE  RATIO 

EXP1=(GAMMA-1 . ) /GAMMA 
EXP2=1 ,/EXPl 

PRCRIT=< 2 . / ( GAMMA+1 . > >*#EXP2 


C - THROAT  VELOCITY  AND  SPEER  OF  SOUND 

VMULI=64.348#EXP2 

VEL=1.-PRA1I0**EXP1 

VEL=SOR  T  ( VMULTKPINT  *144.  *VEL/DENS  TV ) 

VGOUND=l . -PRCRIT  **EXPl 

VSOUND=SORI  (VMIJL  I  *P1NT  *144  .  #VSOUND/DENS IY  ) 

VMAGH=V£L/VSOUND 

C*******#******#t#*t*###****#******#****##******************t**********<I 
C****l  INI  UAL  FLOW  RATE  CONDITIONS  ***«C 

O - VALVE  ONE  FLOW  RATE  CALCULA f ION 

IF  (VAl  VS1.EO.O.  >  GO  10  205 
IF  (DEL  I  A. LI . CRACK 1 )  GO  TO  203 
IF  (DEI ( A . G I . CHANG1 )  GO  10  200 
A5- A  l 
AA-A2 
GO  U)  201 

200  A5-A3 
A6-A4 

201  RA  IEVl=EXF'(  A5I  A6*AL0G(DELT  A  > ) 

ORA I E1=RA I EV1/60 . *VALVSt 

RA IEV1 -ORA  IE1*D£NSTY 

C - VALVE  I  WO  FLOW  RATE  CALCULA I  ION 

203  IF  (VALOS2.EO.O. )  GO  10  210 

IF  (DELTA. LI . CRACN2)  GO  TO  210 

IF  (DELIA. GT .CHANG2)  GO  10  206 

D5-D1 

B6-H2 

GO  10  207 

206  B3=B3 
D6-D4 

207  R A IEV2=EXP t B5+B6# ALOG ( DEL I A ) ) 

RA  I  EV2=RA  IE02*MASS/6(> .  /VOLUME* VAL  VS2 
OKA IE2-RA I EV2/DENS I V 
210  KA I b  U~RA IEV1 1  RA I EV2 

C- - FtuER  FLOW  RAIL  CALCUt  A I  TON 

IF  (Ffl  IRS.EO.O. )  GO  10  220 
KA  I'EF=(‘l +C2*DELT  A+C3*I*EI.  I  A**2+C4*IlEl  I  A**3 
RA  II  F  -RA TEF *AF  1LT  *F I L T RS 
RA  I EF  “RA  I  EF*MASS/6(> .  /VOLUME 
ORA  rEF  -KATEF/DENS I Y 
L-  ---  ADD  IN  LEAK  CONTRIDII I  ION  IF  ANY 
270  IF  ( ALEAK, EU. 0. )  GO  10  230 

RA  IEL=DCOEFF*ALEAK#SOR  I  (  64 . 348*MASS*DEL  IA*1 44  .  /VOL  UME ) 

ORA  I FL -RAT EL/DENS  I Y 
730  RA ! E -KA I  EOT  RA I El  I  RATE F 

--Olirpl.il  INI  IIAL  CONDI  I  IONS 
100  WR I  It  ( NOU  I  r  92 1  ?  >  GAMMA. PRCRT  I  *  VSOIJNR 
WRI I E  ( NOIJ 1 . 9215) 

WR1 IE (NOU T .923 6)  TIME. PEXI rPlNI .DEI  IA.MASS.RAIE 

WRI TEv NOD  12. 9221)  I IME.PEXT.PINI . PRAT  10. DEL TA .VMACH. DENS I  > . 

*  ORA  I  E  J.RA  I  EVt  .LIRA  I  E2.RA  I EV2 .  ORAT  EF  .  RA  I  EF. 

*  ORA  I  EL  .  RA  TEL 

COttttt***************************************************************^ 
Ctput  INTEGRATION  OF  DIFFEREN I l AL  EOUATION  *#***C 

Ct*t* I************#********#**********#********#*****#***************#**' 

C - INITIALIZE  VARIABLES 

H2- I STEP/2. 

PCNT=PINT*144. 

n-=o. 

12=0. 
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i 

2 


T3=0. 

T4=0. 

C - INTEGRATION  LOOP 

500  PRESS-=PINT 
1 1MES-  I  [ME 
KOUNT  =0 . 

C - INTEGRATION  POINT  COUNTER 

540  KOUNT  =KOUN  1 1 1 

T - INTERPOLATE  TABLE  FOR  EXTERNAL  PRESSURE 

550  TF  (TIMES.LE.PT 1ME( INDEX) )  GO  TO  560 
INDEX=INREX4  1 
GO  HI  550 

560  'F  t  I  TMES.GE.PI IML(TNDLX-J >  )  GO  TO  5/0 
CNDEX=INDEX-1 
GO  TO  560 

5/0  !  F  (  INDEX. GT .MAGNI ! >  GO  TO  1000 

FRACN"( T  IMES-P  TIME  <  INIiEX-1  '  )  '  t  P  1  1  ME  <  LNDfcX '-P!  I  ME<  I NOEX-1 '  ' 
T'F  X  !  =F'REX  I  ( INDEX-?  )  I  (PREXI  ( INDEX  )  -PREX !  (  INDEX -1  )  )  YFRACN 

V - CAITUlATE  THE  INTERNA!  GAS  MASS  &  URATE  CORRECTION  FACTOR 

MASS--PRESS/MUL  T 
GENS1 Y=MASG/VOl  UMI 
VCUKR-SOR T ( PA  I  M/PRESS  > 

ocuRR=sem < ucorr  > 

c--  -T  ALOUl.ATF  HE  PRESSURE  DIFFERENCE 
1'EY=0 

UEL  I  A=(PRESS-F-EXI  >/144, 

IF  (I'EL  IA.LE.O.  )  KE V  =  t 

fl-  (  DELTA.  EE.  0.  >  DEL  I  A-0 . 000000 

L - lALCUIAIE  PRESSURE  RATIO  AND  THROAT  VELOCITY 

PRAT  IQ^F'EX  T/PRESS 

VF.L-1  .-PRATUmEXPl 

IF  t  VEL . L I . 0 . )  VEL  ~0 . 000000 

VEL  -SORT (OMUL  T/PRESS# OEL /HENS I Y ) 

IF  (VEL.G! .VS0UN1D  V£L=VSOUND 
VMACH=VEL/VSOUND 

C=-==r-=VALVE  ONE  CALCULATIONS 

C - CHECK  WHETHER  VALUE  ONE  IS  OPEN  OR  CLOSED 

IF  (VALVS1.EU.0. )  GO  10  580 
IF  ( DEL TA . L I . CRACK 1 )  GO  TO  580 

C - VALVE  ONE  IS  OPEN  -  ROW  RAIE  CALCULATION 

C - CHECK  FOR  CHOKED  FLOW  CONDITION 

IF  (VMACH.GF. 1. )  GO  TO  57? 

IF  (DELTA. GT .CHANG 1)  GO  TO  575 

A5=A1 

A6=A2 

GO  TO  576 

575  A5=A3 
A6-A4 

576  ORA TE1=EXP(A5+A6*AL0GC DELTA) ) 

ORA  TE1=0RAT  El/60 . YVALVST  #GCORR 

C - CHOKED  VALVE  -  FLOW  RATE  CALCULATION 

57?  RATEV1=QRAIF I*DENST Y 
GO  TO  585 

C - VALVE  ONE  IS  CLOSED 

580  RATEV1-0.0 

C~~  -==VALVE  TWO  CALCULATIONS 

C - CHECK  WHETHER  VALVE.  TWO  IS  CLOSED  OR  OPF.N 

505  IF  (VALVS2.EU.0. )  GO  TO  5?0 

IF  (  DEL  T  A  . 1. 1  , CRACK2 >  GO  TO  5?0 
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C . -VALVE  TWO  IS  OPEN  -  FLOW  KATE  CALCULATION 

C - CHECK  FOP  CHOKED  FLOW  CONDI  F  ION 

IF  < VMACH.GE. 1 . >  GO  TO  509 

IF  ( DEL  F A . GT . CHANG? )  00  FO  5 86 

R5=R  L 

B6=B2 

00  70  587 

586  B5=B3 
B6=B4 

587  ORA  I  £2=EXP( B5+B6*AL0G  ( DELI  A  > ) 

QRA I E2=QRA TE2/60 . *VALVS2*UC0RR 

C - CHOKED  VALVE  -  FLOW  RATE  CALCULATION 

589  RAIE02=(7RAr£2*DEN8-  v 
00  TO  595 

C - VALVE  TWO  IS  CLOSEI 

590  RATEV2=0.0 

C=====FILTER  FLOW  RATE  CALCULATION 

C - CHECK  FOR  CHOKED  FLOW 

595  IF  (FILTRS.EO.O. )  GO  TO  597 
IF  (VMACH.GE. 1 . )  GO  TO  596 
GRArEF=Cl+C2#'jELTA+C3*DELrA**2+C4#DEL1  A**3 
IF  ( KEY .  EC? .  i  .  >  ORA  TEF=0 . 00 
C?RATEF=ORATEF*AFILI*FIL  IRS/60. *QCQRK 

C - CHOKED  FILTER  -  FLOW  RATE  CALCULATION 

596  RA  TEF=CIRATEF#DENSTY 

C,==i--=0RIFICE  FLOW  RATE  CALCULATIONS  -  ADD  IN  LEAK  CONTRIBIJ I  ION  IF  ANY 

C - CHECK  FOR  CHOKED  FLOW 

597  IF  <  ALEAK.  EO.  0.  )  GO  10  .599 
IF  (OMACH.GE.l.)  GO  10  5*>8 

RA  1  EL=DCOEFF FALEAK*SORT  ( <s4 . 348*MASS*D£L  I  A*1 44 .  /UOl-UMF  > 

ORA TEL=RATEL/ DENS  1 Y 
GU  (0  599 

L - CHOKED  CONDITION 

V'H  RA!EL=(?RAfEL*D£NSl  Y 

C-- CALCULATE  THE  TOTAL  FLOW  RATE  AND  PRESSURE  DERIVATIVE 

59°  RA  1 E =RA  IEV 1  ^  RATEV.'M  RA !  EL  IRA  1 1 i  F 

F’REDER- -MULT *RA1  E 

l. - PERFORM  INI  EGRA  ■  ION  CALCUI  AT I  ONS 

600  GO  10  C 650  t  700 » WOO t  900  >  *  KOUN I 

650  l  l-ISIEPJTPREDER 
PRESS- PIN  NT  t  '2. 

T  IMES- IIMEFH2 
GO  TO  540 

700  12= i S TEP*PREDER 

PRESS=PTNT  +T2/2 . 

1  TMES-I IME+II2 
GO  TO  5A<) 

8'*0  1 3“  I  SI  EP#PREDER 

PRE8S=PINT+T3 
T1MES-I IMEI (STEP 
GO  FO  540 

900  T  4= I ST  EP*PREDER 

C - CALCULATE  NEW  INTERNAL  PRESSURE 

K0UNT=0 

PIN  r=F’INI  +  (  III  2.  *T2+2 . *  T3I 7  4  > /6 . 

C - CONTINUE  THE  INTEGRA  I  I ON 

1000  1 1 ME=TIME+ I S IEP 

IF  < I  IMF .GE. ! IMEND'  GO  10  1500 
L INDEX=L INDEX+1 
F;RMARK=L  INDEX /DIVIDE 


*  * 


1500 


LMARK-RRMARK 
RMARK -LMARl\ 

IF  ( RKMARK . NE . RMARK )  GO  10  500 
WRI 1 E  <  7  r  V21 A )  fIME 
PINT=F'INT/144. 

PEXT-PF.X  1/144. 

WR1  IE (MOOT  rVE’A)  TIME » PEX  !' r  P  IN1 1  KELT  A  r  MASS  -RATE 
URI  IE ( N01IT2 r  9221 )  TIMErPEXT  »PINT rPRA'1 TO  j  DELTA? VriACH* DENS  TY  r 
ORATE  1  ? RATE01 1 ORA  TE'.' *  RA I  FV2? HRATEF .  RAT EF  ? 
lIRATELrRAIEL 

P1N1=PINT*144. 

IF  (TIME. I  I.riMENLi)  Ou  10  50Q 
STOP  'END  OF  IN ! EfiRAT  TON' 

F  Nil 
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Initial  flow 
conditions 
(Q's  &  m's) 


300 


Output 

initial 

conditions 


INTEGRATION  LOOP 
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Kount=1 


650  | 

i  2 

700  t.  _ 

Calculate  T1 
set  up  for 

T2 

i - , - 

Calculate  T2 
set  up  for 

T3 

INTERMEDIATE  STEPS  IN 
INTEGRATION  SCHEME 


- - PROGRAM  E'RESSM .  FOR 

- THIS  l!.‘  AN  ATTEMPT  10  PREDICT  THE  INTERNAL  PRESSURE  UE  A 

. -SOUND  TNG  ROCKET  AS  IT  ASCENDS  THROUGH  [HE  ATMOSPHERE. 

- THE  PROGRAM  CALCULATES  THE  INTERNAL  PRESSURE  HIS  I  (TRY  OF  A 

—  —  SF'T  OF  MULTIPLE  INTERCONNECTED  CHAMBERS .  AT  PRESENT  THE  PROGRAM 

■ - AltOWG  ONLY  TWO  VOLUMES,  THE  PROGRAM  INTEGRATES  THE  DtFETRENt  I  A! 

- EOUAIIONO  NUMERICALLY  USING  A  FOURIH-ORHtR  RUNGE-KU T1 A  SCHEME. 

-—-THE  INTEGRATION  IS  HONK  IN  THE  " REGS’  HUDROU'JNE  CALLEIi  BY  THIS 
--  -PROGRAM.  IMF  UTFEERENIIAL  E'OUA  i  IONS  ARE  EVAIUATT'D  AT  OAR  LOUS 
---I  IMIS  BY  THE  ’PRIONS’  SUBROUI INF  WHICH  IS  LAI l FP  BY  ’PROS.” 
-author:  *■',  KREBS 

0  I  Ml  NS  1 1  IN  AUX  <  S  r  P )  r  I  I  M  I  ISC/'.  PRESS  (  '  7 1  ‘RE  HER  (  P  ) 

TUMI  NS  I  ON  I  Ab(S>  7  I  L  T  LI  'PC,' 

PIIMCHIN  -  I  Ml  M)  IS/  OAT  OS  t . CRACK 1  ■'  HANG  I  .  A  !  -  A.  .  A  1  •  A  4  . 

■Y  OAT  VOP.CKACKP-CHANGPrHl  ,  B.’  r  BA  *  B4  7 

Y  OAI  OS  < » CRACK  T - ‘'ITANGA  ■ 1 . 1  -r:M\?*C4» 

Y  OAI  OS4.T  KA<  K-t.l  HANG4  7  H  1  -BP.TU  014. 

Y  !  1 1  I S  l  .  AI  H  I  1  r  F  I  r  E  ’  7  E  A  -  E  4  .  A'  K"AK  I  t 

*  !  !  I  I  SP .  Ai  1 1  T  P  r  f  1  -  F  7 1  1  - 1  4  .  *V  E  AK  > 

i  OMMIJN  OUTPUT'  I  I  Ml  -  P  !  N )  1  •  T  E  X I  r  -  RAT  l  .BE  l  *  AI  .HI  Nb  l  -MASS  I  ?OLI  1  s 

*  OMAt„-l!  -RAII-  I.Rh.sO!  -  RA  »  E  Op  .  RA  Iff  I  .  RA  T  E  LI .  ORA  I F  I 

Y  DKA'E.MIKAIE)  rliRAtl  I  .  F'l  N I  ;> ,  ORA  T .  PRA  I  T  i  DEI  TAP? 

*  DLL  l  A,i»HF.NSPr  MAGS1*  .01  L  P  -  VI  L  -5  7  VMAT.HO  r  OMaCFLT  . 

Y  RA  I E  .»  t  KA  ( l  5 ,  RA  i  F  OS  .F'A  I E04  .  RA  I  E  F  J  r  RA  l  E I  -  OKA  T  M  « 

*  ORA  l  E  4  .  ORA  II  '*.  ORA  1 1  I’ 

l  OMMON  /REEF  HO  '  001  1  1 00’. .  MU!  F1.MUI  I  ’-VEIM  T  -  EXP  l  *  OS()l  INC  • 

*  PRE-XT  (  1 00  l  <pl  IKK  1  AO  >  -  INBE  v-iK.Ot  E  T  » MtfON1  ! 

REA'  !  IMI  ISrMASSI  n-ASS  'rMUl  11.001  !  .' 

EXTIRNAi  f ’NEONS 

t'A  I  A  ORA  I T  t  r  ORA  I F  MIKA  1  t.T.WRA  T  L  4 .  ORA  ! 1  I  .t!PA  f  I Pr  ORA  H.  1  -  ORA  JL  2  r 

*  RATE  [.RATI  .> .  RA  I E  ,<  r  RA  H  F  L. RATIO  PrRAll  »  '  rRAIEl  .’rRATFVI  , 

Y  RATI  Ol’.RATF  OX  rRA  I E  04  .  BCOE  M  .INDEX. I  I  NUF  X  .  LMARK/1  9*0  .  .  0  .  VO  • 

Y  .’.'BIT/ 

-  -FORMAT  SI  AT  l  MINIS 
9000  FORMAT <P0a4T 
9100  FORMAT l 1  I  Or  7f 10. 0) 

9 10  I  FORMAT <01  10.0) 

Rpoo  FORMAT (/lOXr 'HAS  PROPFRT  TEST '  / 

*  t  f»X  r  '  I  YE'E. !  *  r  SA4 

*  l  OX r ’ MAIN  OUI  UML  =  'rMO.P.'  CU  FT'/ 

*  I SX. 'SECOND  VOLUME  -  '-ItO.I’r'  CU  ET'/ 

Y  IbXr'INTICAl  PRESSURE:  -  'rElO,'.’.'  E’SJ'/ 

Y  I  S  X  »  '  IF  MPERA  1 URF  -  '  .  T  1  0  .  I  .  '  DE  GREF  S  E  '/ 

t  1SX. 'GAS  CONSTANT  '-EIO.I’t'  I  T -I  R 'I  B-Dh C  R ' 

Y  V) 

9. *01  E (IRMA )  ( /  LOX r  '  MAI  VE  ONE  PROPER!  I  I.S ;  '  / 

Y  15X . ' T  YPE l  ‘  . P0A4  F 


r  SA4 '' 
r  I  10..’. 

■  i  lo.;* r 

’  .  E  1  0  ,  ? . 
.110.1. 
-  e  i  o .  ;>  j 


CU  FT'/ 

CU  El'/ 

E'SIV 

DEGREE  S  E  '/ 

I  T--I  B'l  B-BEC  R' 


9 .’OP  FORMAT  <  1  b.x  .  'NUMBER  OF  RE  l  TEF  VAl  OLS 

Y  I  bx ,  CRACK  I NG  E’RESSTJRF 

y  j sx . 'i.ukvi  change  point 

Y  i:./7  'l.OE.I  I  ICLE  NT  1 

Y  TbXr'COF.FE  IC1E  NT  .* 

Y  I SX . ' COt E  E I  0 ItN T  A 

Y  tOx  7 ' CUFEE  CLIT  Ni  4 

9P0P  FORMAT t.'lOX 7 'VAl  OE  TWO  PROPERTIES:'/ 

Y  iSX.'IYE'i: 

9204  EURMAli  '/IOX.  'VAl  VE  (HKIF.  E'KOF'F.K  TIES: 


* POA 4 > 

-E 10,0/ 

' ,  i  to,  ;> , 
.EIO.,'.' 
-  E  I  0 . 9  / 

,  E  1  O  *  A  / 
.E  I  (; .  X/ 
'EV',,T/  ' 

7  .’0A4 ) 


PS  I  '  ' 
PS  I  '  / 


9 .’OS  E  OkMA  I  l  /  I  Ox  *  '  VAl  Vtr  FOUR  PROPER  !  Tf  G : 


■•’.’06  I  IlKMAI  (  '  I0X. 'F  II  H  R  ONE  PROPER I  It-'!!  / 

*  LSX.'IYPL:  '  / 20A4  > 

920/  FORMA) <  i;.X»  'NUMHEK  OK  F  II  I t  Kh  -  yMO.O/ 

*  IbX.'EXll  AREA  --  ',I10,S.'  SO  IN '  / 

*  1 1">X  i  '  CUFFF  TCIFN  F  1  -  'rFJO.X' 

*  IOXt  00EFF1C1ENI  2  -  'riw,.i/ 

*  lOXr  'COfcn  ICltNI  X  -  ',ELO..V 

*  10X,  '('Otf-KICtH'ir  4  ',MO.X/> 

9200  !  OKMA  I  (  /  l OX »  ' HI.  It R  I  WO  ■•'ROPKR  ( J  F.S  !  ‘  / 

*  1SX .  '  I  YF’F !  '  *  20A4 ' 

9 'JO  9  FORMA!  </tOX,  'HOUR  I  EAR  ONE  PROPF  K  I  ft  H  S  *  / 

*  IOXt  '  I YPE  (It  SEAL:  ',20A4> 

9210  KOKMAK  IHXr 'EFFECHOF.  AREA  -  '.110.8,'  SO  Fl'/> 

99  1 1  FORMaI  (/ IOX,  'IKIOK-  l  EAK  f  WO  PROPER  I  l£S ! ' / 

Y  LOX-'IYPt  (It  >U  AL  :  '  * 20A4  > 

9939  FORMA!  <//l()X. 'CHOMNH  PROPER  1  ft.S !'  / 

t  ISX, 'RATIO  Ot  SPEC  frit,'  HF.A  I S  -  M  10.1/ 

Y  lf»Xr  'OKI  f  ICAL  Pkt-SSORt  RAffO  --  'rFt0.4,' 

Y  tOXr  SPtt'fi  Ot  SOIJNH  -  '.FlO.i.'  FES  ///' 

'•>  >1  <  I  OKMA!  (  MAIM  001  OMt  (,Al  COLA  I  IONS'  /  ' 

99  I  4  f  OKMA  I  <  ’  SFUINHAKY  VOLUME.  CAECUL  A  I  f  tlNS '  /  ) 

99  I  0  K  OKMA  I  I  XX  r  fXIKRNAl  I  N  1 1 RNAI  ' , I 4X* '  INI  t.RNAl  '  -  4X  r  HUAI.  MASS'/ 

*  !,X  .  I  IMI  ORtSSURE  '  ),XX,  •  IlH  FA  P'-'lX,'HftS  MASS'. 

Y  l.Xy'tlOW  RAH  '/SXr  'SECS'  yAX,  PS  l  '  .  9  (  /X  .  '  PHI '  '  r  9  <  ,  '  I  KM '  , 

Y  OX.  IKM/SFT'/) 

9  J*>  I  ORMAI  H-  9.  1  .2  !  0 . 2 .  f  1  0 .  ,  F  l  X . !  > » F 1  4  .  /  > 

9  ’  I  '  I-  (IRMA  I  (  I  V  X  »  P  R  t.  S  S  (I  K  t'  rllXy'PRI  SSUKt  0  H  t  '  .  OX .  '  1 01  AL  MASS'/ 

*  SX ,  '  I  I Mf  ' . AX .  tXf  MAIN  OOL  SH‘  001  ' -  OX , ' MAIN ' , 

*  .<>V.  fX;  '.AX, 'll  OW  RAfl.'/SX.  SICS' , AX, 'EM ' ,4</X» 'PHI '> ■ 

t  /X,  I  PH/St  O' /> 

9'MH  f  ORMAI  <1  9.  I  .  ¥  U). 9,93  10 .  ,5 , 1 1 4  .  / ' 

999 A  f  ORMAI  (SAX,  'OALOL  (INF '  ,  9X , '  0A1  Of  I  WO  .OX.  FI!  I  KR  ONE" '  » 

*  9Xr  t.t  AK  ONE.' /t,5X.  'P  R  t.  S  S  0  R  £  '  ,  1  ^ .  4  ( 9X*  '  tl.OW  RAIL1'/ 

Y  ">X,  I  fMt.  f-Xl  INI  RAIMI  UlFIR'rlXr  MAOH  OtNST  I Y  ' 

*  4  (  '  001  OMt  MAS!.  /2X,  '  SECS'  */<  IX» 'PS  l '  )  ,  ,3X .  '  PE/PI  '  , 

*  PS1  NO  I  KM/CD  FI  '  ,  4  (  '  OF/S  l  CM /SI- C  ')/> 

999  I  FORMA  I  f  F  A  .  I  ,  I  X .  2f  A .  2 « F  / .  4 . F  V .  X , h H .  X , F 9 .  S. 4  <F W .  4  . F 10 .  A  >  > 

9999  I  ORMAI  <9X? 'MAIN  001.  REE FERENCE '  , 4\ ,  E  X  1  ERNAi  REFERENCE  , 

*  I4X.OAIOE  IHRtF' . /X, 'OAl  OF  FOUR'  ,0X,  'F  fl.  I  £K  IWO'r 

Y  9X, 'I  EAR  IWO'/tlX, 'PRE SSURF ' - 14X, 'PRESSURE' ,15X, 

*  4  ( 9X »  now  RAlt  > /2X  » '  I  TME  '  ,  1’ (  XX ,  '  RA  I  10  I .  E  t  F  ft  MACH'), 

Y  .’X,  '  Of  NS1  I Y  '  , X (  '  001  OMF  MASH  '  )  ,  ' OOL.IIMF.  MASS'/'XX 

Y  'SECS  PI/P2  PS  i  NO  PE/I  -9.  PS  l  NO  IJiM/UJ '  , 

*  '  F-  »  '  .  4  <  '  Of /S  LPM/SFC  ')/) 

922X  f  ORMAI  (FA. 3  , 9  < FE! .  4 , 9S7  ■  3  >  ,  F  9 .  Sr  4  ( FH .  4  ,  F  10 . 6  >  ) 

- .  iMt  KEF  ERF  NOES 

i  I  ME -=0.00 

C - -At  I.OCArt"  OAl  A  FfLfS 

N1N-X0 
NOIII  =NINI 1 
N0IJI2=NUUr+l 
Noiirx=NUoi  f:» 

MOil  1 4-NOll  I  IX 

C  . LNPU1  ANH  OOIPO!  OF  OflAL  INFORMA FION 

(• .  CAS  PROPER!  fES 

KE.AIK  N  1 N,  °000 )  1  1  1 1  F 
PE  All  (N1N,  9000'  HAH 

HEAP  l  N  f  N ,  9101 )  OOL  t  .  VUI.2 ,  r'TN  I  .  I  EMP .  RHAS ,  HAMM  A 

WK,r.b  (Nour./ooo)  ixilf 

WR  TIE  (NOIII  ,9200)  HAS, 001.1 .001  2.  PIN  F  .  I  FMP  ,RUAS 


UK'  1 !  E  <  NO  l!  12,9000)  I  It L E 
wk >  1 1-  <noui2.9213' 

UP  l  S  E  ( NOU 1 2 , 9220 ) 

it  (voLi’.tti.o.  >  on  ro  5o 

UK  1 1 E t NOU T3 . 9000 >  r  I  n  K 
UR1  ffc  (NOUM.VOOO)  TIM  P 
URfTP  (NOlit  4.9214) 
wrt  it  (noth  4.  y 

-h.VP  001  (IMP  PkEftftUkt  SYftlFM 

-  -<>Am  ONfc  PROP'EkftES 

P<  A0(N1N.9000>  :  U!  t 

K'i  AlUNIN-9101  >  OAI  OK  L .CRACK  1 , CHANG  I  rfl  i  .02.03.04 

i  rt  '.NoiH  .yyoi  >  mi  t 

!•  I  OAI  00  l  .tll.Oi  )  00  It)  10O 

I  f  I  E  ( NOD  1  .  '■>209 )  OAI  OKI.1  .'HACK  1  -  CHAN  ft  J  .At  ■  AI-i .  A3 .  A4 

-  -  "tii  oi  run  pkop-fc k i  its 

1  RP  Am  NIN.  9000'  1IILE 

k:  AIKNIN.ytOt )  OAI  VS2. CRACK.’. CHANG2,-'H -H9.H3.H4 

Wkl  I  E  (NOII I  .9903  1  I  HIP 

II-  <  VAL009 .  Ell  .  0 ,  >  00  10  110 

Uk  f  rp  ■.  NOU  I  ,  9202 '  OAI  0b2 ,  CKACK2  . 0HANG9  .Ki.R2.H3.  04 

. i  ILTER  ONE  PROPER  1 J  Eft 

1  1 0  RE  Alt  ININ.  VOOO )  I  l  PLE 

RE  AH <NTN» 9101 )  F1LTSI.AFIL 1 1 .El .E2.E3.E4 

UK' ME  (NOU  I  .9204)  imp 

IP  <F1L ISl.EO.O. >  GO  (0  120 

Wkl  rE'INOIII  .9207)  Hi  TSt.AFILri.EI  .E2.E3-E4 

- IE  AN  ONE  PROPER  TIES 

190  REAM  NIN.  9000)  El  IKE 

ET-AIt  ( N  IN .  910 1  )  ALP  AM 
UR1IF (N0U1 .920y>  I  HIP 
IF  (AIEAM.PO.O.)  GO  10  130 
Uk  J  I L-  ( NOU  I  .  921  0  )  At  P  AN  I 
- -- sFCONn  OftLUMF  PRESSURE  SYSIEM 
I  30  [>-  (00L9.FU.0.)  (.0  TO  170 

- OAI OP  T  HREE  PRUPF  R I TES 

RE  AlHNTN.yOOO'  TITLE 

RPATKN1N.9101 >  OAT  0S3.CRACK3.CHANft3.Ct .(  9.C3.C4 
Wkl  IL'(N0U!3.y204)  I  [  fl.f 
IF  (0AL0S3.FQ.0. >  GO  TO  140 

Wkl 1 M NOU T3 . 9202  >  OAI  OS 3 . CRACK3 - CHANG  3 . C 1 . c2 . C3 . C4 

- OAl  OP  POUR  PROPP k T i Eft 

1  40  Rt AU<N IN. 9000  >  llllf 

RE All  < WIN . 9 1 01 )  OAL  0S4 . CKACK4 . CHANG4 .01.02. 03 . 04 
Wkl  IF  ( Null  1 3 . 920ft )  UUP 
1 1-  <  0AL0S4  ,  Ell ,  0 .  )  GO  TO  1 50 

UR  r  I  £  (NOU  I  3.9202 )  OAI.  0S4 . CRACK-* . CHANG4 .01  .02.05.04 

-  --Ml.  I P  R  1UU  PROPER! IE S 

I  50  KP  AO <  NIN . 9000 )  I  TILE 

REAM PIN. 9101  )  I-  U  l  S2 .  AF II.  12.1  1  .F  9.1  3-P  4 
WR>  T  E  (  NOU  I  3  ?  9900 '  f  T),l  E 
I  (F TLTS2.E0.0. >  GO  10  1  AO 

WR  l  ( P  <  N(JO  T'3 . 920/ )  F  1 L  CS2 .  AF  U  !  9 .  F  1 . 1-  2  -  K3  -  P  4 

• - LPAK  TWO  PROPERTIES 

160  RFATuNtN.9000)  TITLE 
RF AO i NIN. 91 01 )  ALFAK2 
Wkl 1E(N0UI 3.921 1)  TITLE 
IF  (ALEAK2.E0.0. )  GO  TO  170 
Wkl I E ( NOO ! 3.9210 )  ALEAK2 
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•INPUT  OK  EXTERNAL  PREf  '  HtSIORY 

REAH<NINr9101 )  TIMEN  1  -^DIVIUE  .LIMIT  S<4> 

REAR  (  N1N  .  9.i\/b  J  MHtii- 

REAIKNIN.  910J  )  (PI.  .  REX  I  (N) .N=l »  MAI3N1 1  > 

CONVERSION  10  PROPER  .S 

I EMP=TEMF’+459 . 67 
PIN !  =PINI +PREXT  ( 1  '144. 

CAl  CULATION  OF  INI  .L  CONHIIONS 
MAIN  VOLUME 
INIERNAL  OAS  LiENSlTY 
M  ■■  11 -R(3AS*lEMP/V0Lt 
MASSl=PtNI*144./MUL!l 
MENS  I =MASS1/V0L1 

PRFSSIIKE  1’IFFERENI  IAI  ANP  PRESSURE  RAItO 
PEX  I  -"FF:EX"i  <1  )/144, 

UU  l A l “PIN T— PEX I 
PRA I  ]  -F'EX  I  /P  (NT 
CRITICAL  PRESSURE  RATIO 
I-  XPl  ~  <  OAMMA- 1  .  1/OAMMA 
FXK>-t./EXPl 

PRCRI'I  -<2.  /(OAMMAl  J  .  >  >**EXP<? 

IHROAI  VFL001IY  ANU  SPEF.U  OF  SOUNIl 
VMUL  I  -A4  .  S4BKF  XPL' 

VL l  I -1  .  -PRA ! 1 **FXPt 

VFl  l-SURHVMO!  T*R  [M  I  *144  .  *VEL1 /lit  NS1  > 

VSiiUNIi^J  ,-PRCRl  l**EXPl 

V'>OIINli=SOR !  (VMUl  I+PJN1Y144. fVSIJUNU/IiENSl > 

"MAl.HI-VEL  1/VSOUNH 
SFt  ONI]  VOLIIMF 

I  F  •.  VOL  J  .  Ell .  0  .  )  00  10  /!>■.) 

MU!  I  J-RfiASt  I EMP/  VOL.’ 

MASOL' -PIN  1*144  .  /hill  II' 

OENO.'-MASSr’/VlU 
REt  !A?=PIN!-PIN! 
tit!  IA.<=-P]NI  -PEX  i 
PRA1  ,'=-plNl  'PIN! 

FRAli-PF.XI/PtNI 
VFL.’-  1  .  -PRAK'+*FXPt 

VFl  t'-SQRI  ( VMIJl  I*RtNI'M44.*VEL'.’/IiENHL,> 
vmach?--ve  LtVVSOUNIi 
-OUTPUT  INITIAL  CONDITIONS 

WRlIEtNOUT  r  9?  1 ?  >  OAMMA.PRCRU  rVSOONIi 
UR  1 1 E  <  NOU  T  r  9213 ) 

UR  l  IE  (NOU!  jyt'lt.) 

UR  f  I E  <  NOt  l !  r  V„>  l  A  )  I  l  ML  .  Pt  X  f  » P I  N  1  .  OEL  I  A  I  -  MASS  I  ,  RA  I E 1 

UR  1  rt  <N0UI?-9.!L1  '  I  IME.PtXI  .PINT  >PRAIlr!iEl  I  l « VMACHt  .  PENH]  . 

'  (TRAIT  1  .RAIEV1  .  OKA  1 1  2  .  RA  T  F  V2  •  ORA  I F  1  -  RA  I EF1 . 

>  (TRAIT  1  .RAT  El  1 

rt  ( V0L2 .  E(T .  0 .  '  bO  10  400 
UK  L  IF  t  NOU  T.(  -y?  I  4  ' 

UR  I  IK(N0UT.T»9.>I  V) 

WR  I  Tf  I  NOU  IS*  91?  1 0 '  I  IMT  -It  x  i  ,P1NI -PINI -UK  lA.’rUEt  !A3»RAIF.T 
URI  IF  (  NOU  I  4  .  9.’?S  >  I  IMt  »-'PA  I 2.  Ill  I  I  A,’.  VMA0H2 .  PRA  1 7,  r  DFl  I  A.?  -  VMACHS  t 
I  DENS? .  ORA  *  E  ?  .  RA  I  E  VS  r  ORA  I  E4  .  RA  I F  V4  -  OKA  If?, 

i  RAIN  ?.ORAII  ’».RAIEI  ? 

bl  I  UP  I  UR  IN  I F  ORA  I  l UN  l  OOP 
-IN!  I 'AT  VAllltS 

PR!  SS(1  >-=PTNJ*F44. 

PR!  SS  ( ? )  -  F'RF'HS  (  t  > 


*  * 


-  X N !  1  OKA  I  TUN  LOOP 
V.>0  l  IHl  IS(  I  )-I  IMF 

Int  IS(2»  =  I  IMF.+  ! b ( fc.P 
I  X  ri  X  I S  <  3  >  =  l S 1 £P 

- RFSh.l  ERROR  Mt’ltiKlii 

PKF  nth- ( 1 '  =0 .  bO 
PKL-HERt2>-0.b0 

- I  AL  l  INTEGRA! TON  ROUNf INF 

I  A!  L  RKGS  (LIMITS.  F’KESS • PKEIltR  i  2  .NSEG  i  .  PRhONS  .  AUX  > 

II  (  11  ME. HE.  rihENri)  GO  10  AO0 
I  XNriF.V-LINLlEXU 

RRMAKK-I  INIiEX/riIUXHt 
!  MAKR-RRMARK 
KMARK^I.  MARK 

IF  <  RRMARK , ME . RMARK )  00  >0  000 
A 00  WR  U  E  7 , 92  J  6  )  i  IMF 

UR X  rElNOIM  *?C16>  I  IMF  rPEXT  .PINT  1.JIELT  AX-MAHS1.RATE5 
UR  l  rF  <  NOU !  2  7  0221  )  X  LME  »E‘EX  I  -  PIN  f  1  r  RRA  1 1  .  UELTA1 OMACFM  ,  IiENbl . 

ORA  I  £  1 r  RA ( F01 . OKA  t  E2 . RA ( F 02  r  ORA !  F 1  .  RA  TEF  l . 

ORA  I L 1 t RA Fft l 
IF  ( 001.2. F 0.0.  )  00  10  /OO 

UK  X  (F  tNUllI  3.0210'  I  lot  rPEXT  rCXNI  I  rR  IN  I  2*  Uhl.  I  A2.  LiF.l.  1 A3 .RAT 
UR  I  1  F  ( NOU I  <1*02 23 )  1  i  ME .  PRA 1 2 .  lit  L  1  A2  r  0MAI.H2  •  PRAT  3  -  IiF.1.7  A3  .  OMACH3  - 
*  IIF.NS2 . ORA  1 E3  •  RATF03  -  ORAT  E4  ,  PA  FE04  .  ORA  !  F  ? . 

i  RATI-  F  2  -  ORAT  I.  ''  r  RAT  EL  2 

'(>o  if  <•  rim  .li  .  i  rMtNio  go  io  ooo 

Slop  '  F  Nil  OF  J  N  TF  ORA  !  >  ON ' 

F  Nl* 
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Start 
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iKtJufihiituiiliiiift  mWWnpfabtmiMittiiniiiii.  Bito*  AnfMihMit  *  t*.*p  ft  nnino'^  hiVi*  »H  i 


Subroutine  RKGS.  FOR 


Modified  RKGS  routine  of  Reference  4;  called 
by  main  program  to  integrate  the  differential 
equations. 


The  following  lines  nave  been  modified  to  eliminate 
the  use  of  the  external  output  routine: 

RKGS  1050 
RKGS  1530 
RKGS  2280 
RKGS  2570 


Integration  Technique:  Self-starting  foui+h-order 

Runge-Kutta  solution  of  a 
system  of  first-order 
ordinary  differential 
equations. 


Subroutine  PRIiQNS.  FOR 

Called  by  RKGS  subroutine  to  evaluate  the  differential 
equations  during  integration  process. 

Major  liquations  Used:  Equations  (12),  (13),  (14), 

(All),  (A  18),  (A 25),  and 
those  of  Table  Al. 
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SUBROUT  I NE  PREQNS  (  T IME I N ,  PRESS  f  PREDER  ) 

- PREONS. FOR 

- THIS  SUBROUTINE  IS  USED  UIIH  THE  'PRESSM*  PROGRAM  TO  PREDICT  THE 

- INTERNAL  PRESSURE  OF  A  SOUNDING  ROCKECT  PAYLOAD.  IT  CALCULATES 

- THE  FLOW  CHARACTERISTICS  OF  IHF.  PRESSURE  RELIEF  VALUES.  FlLIERS  AND 

- ORIFICES  USING  (HE  EMPIRICALLY  DERIVED  HOW  CURVES.  I!  EVALUATES 

- 1  HE  DIFFERENTIAL  EQUATIONS  FOR  THE  INTEGRA! ION  SUBROUTINE  "RKGS* 

--  WHICH  IS  CALLED  FROM  THE  MAIN  PROGRAM. 

- A'IIHOK;  C.  P.  KREBS 

D  !  MENS  ION  F*PE  SS i 2)  .F'REDER*'  2 ) 

COMMON  '  INPUTS/  VALVSI  r  CRACK  1  r  CHANG  1  -  At  .A2-A3-A4. 

♦  VALVS2  ,  CRACK.;  r  CHANG? . B1 « B2  >  B3  r  B4  * 

1  VALVS3 .CRACK3 < CHANG 3  r  1.1  •  l  ?•:  3.C4. 

t  VAI  Vb*'l.fiACK4. CHANG-?. D?  tD2tD5.D4. 

♦  F  1LTS1  -AF  IL  T  l  rF  I  -E2-F.  T-r'4.  ALEAK  1 . 

♦  F  !L  rs:>-  AE  J>  T?.F  '  •  t  .  E  5-E4  r  ALEAK 2 

i  UnMUR  •  y  ■  rr'!l  I  '  I  IMF  .PIN  I  !  rfEXT  ,PRfl  •  1  .DEI  f.Al.DhNSI.MASSl.VELl. 

»  ijr.fi  -  RA  f  *■"  I  •  RA I  £  <'l .  F:A  *  EV? .  RA  TEF 1  rRAIEL  1  jQRATE’  . 

f  UPf*  Tf  2-ftKATF  I  rllRATLI IN ! 2.PRAT 2.PKAT  5.firiLIA2> 

+  DEL  1  A3.  DFNS.’r  MASS.’ .  VEl  2.VH  3  r  MMACH2  .  VHACR3  . 

f  E*A  I E  2  -  RA  I  E  .3 » RA  I F  V  I  ?  RA  !  E  V*1 » RA  TFF2  .  KlA  !  E I  2 .  URA  TE3  * 

t  fl'il. {IRA  Tl  . GRA  1 1  2 

•■(jriMl'N  -'F:FF  t SS/  VULl -V'JI  .’.MUi  !  I  .MUL  :  J.VMUl  !  >EXcl  rVSOUND. 
f  ERL  X  T  <  >  00 )  t  P II MF  •  1 00 '  .  I  Nl'F  X  .  Di  OEEF  .  MAGN 1 ! 

PEA1  Mr- -SI  •  MASS/rMIJL  f  1  -  Mill  1  2 
-  !  RANGE  I-  R  I  HF  f  i  MF 
( I  ME  :  IMF  IN 
■  1  ME  b  -  II  MF 

--  -IRANSFLR  IFF  INIERNAl  PRESSURES 
PRESS!  --’F.ESS I !  ' 

PF:ESb2-c'RESS  ( ' 

PIN!  1 -E’RESSt  1  '  '  1 44  . 

PIN!  2“F'F:ESS ( 2  > / 1 44  . 

- 1  N  I  ERF  OLA  IF  TABLE  FUR  FXTFRNAI.  PRESSURE 

LOO  IF  C 1MES.LE.PT IME< INDEX' )  GO  TO  110 
INDEX=INDEX+1 
GO  10  100 

HO  IE  (TJMLS.GE.rM  IMEtlHPEX-l'l  GO  !U  I  .To 
i  NpF.  < -INDEX—  L 
GO  10  110 

I"'*  !  *•  (  INDFv.G1  .MAgN!  1  ■  SIUP  INDEX  100  LARGE 

FROCN  < 1  !  MFS-P I  l ME  <  :  NfiEX-’  >  '  -  ipl  IMF  <  INDEX  ) -E  I  IMEl  INDEX- 1  >  } 
it  vr^c-PEXf  (  INDFX-l  *  F  *.  PRf  XI  (  1NDE  X  '-C'RFX  !  (  TNDFX-1 )  '  *FRACN 

otmmmi  t  r  rtf  t  t+toitf»M*t$*t***tt*+**t*?*t**t******t**c 

HIT  INNER  VOLUME  CALCULATIONS 

umn  uiiHmutTmTmmmi  tfSTittmTmmsmttTtmmmTC 

11  !  VOL? .  Eli .  0  •  ■  GO  *  O  SO** 

*,l  I'pl  A1  F  GAS  PROPER  I  HS 
MASS? - PRF SS? /MU!  ! ? 

PENS.  -  MYSG-T/VOi 

--  »  cl  .  *  *L  A I  E  THE  ORATE  s  UPpEi  I /UN  r  .H  f  Of. 

< *  !  iRR?'  SOR 1  PRE  X  i  (  1  >  .'PRF  bS? ' 

>.*(  *n--P2=SflPl  ■'  OCORR?  • 

- - !  A:  I1'  ATE  i  ML  PKESSUE'E  f'i^EI  REK'CES 

-  !  F  i“( 

;«  !  I  A.*  •.  PRESSL'-PPtSSI  '*  *144. 

!■!  :  '  A3=  (  PPfc  SS3-I  E  x  1  )  / :  4  /•  - 

I  ‘  >  DE L ! A3 . 1  I  . 0 .  >  Kt  rJ- 1 

I'  '  DEI  JA2.1  E.O.  >  DE?  I A  1-0.000000 
IE  *.DFLrA3.LE:.0.  >  DH  IA3=0. 000000 


J  . 


c - 

240 

c=-=== 

500 


C - 

540 

400 

C====: 


-CALCULA I E  PRESSURE  RATIOS  AND  THEOA!  UELOCI'IES 
PKAT2=PRESS1/PRESS2 
&RAf3=PEXT/PRESS2 

ot-j  2=i.-prat2*#expi 

IF  < VEL2 . LE . 0 . )  VEL2=0. 000000 
VRL2=SGRT  < VMULT#PR£SS2*VEL2/DENS2 ' 

VET.  3=1 .  -PRAT  3#*£XP1 
IT  <VEL3.LE.O.)  VEL3=0. 000000 
V£L3=SGRT ( VMULT*PRESS2*VEL3/DENS2 ) 

IF  (VEL2.GT .OSOUND)  VEL2=VS0UND 
IF  <  0EL3 ,  GT . V50UND >  V£L3=VS0UND 
VMACH2=V£L2/VS0UND 
VMACH3=VEL3/VS0UNB 
=VALVE  THREE  FLOU  CALCULATIONS 
-CHECK  WHETHER  UALUE  THREE  IS  OPEN  OR  CLOSED 
IF  '.VALVS3.EQ.0.  >  GO  TO  240 
IF  <  DELTA2.LT . CRACK 3 )  GO  TO  240 
-VALUE  THREE  IS  OPEN  -  FLOU  RATE  CALCULATION 
■CHECK  FOR  CHOKED  FLOU  CONDITION 
IF  (VHACH2.GE.1. >  GO  TO  230 
IF  (DELTA2.GT . CHANG3)  GO  10  210 
C5=C1 
C6=C2 
GO  TO  220 
C5=C3 
C6=C4 

0RATE3=EXP  (  C5+C.6* hLOG  <  DEL  I A2  >  * 
QRATE3=QRATE3/60.*VALVS3*0C0RR2 
-CHOKED  VALVE  -  FLOU  RATE  CALCULATION 
RA  TEV3=QRA TE3»DENS2 
GO  TO  300 

-VALVE  1HREE  IS  CLOSED 
RATEV3=0.0 

-VALVE  FOUR  FLOU  CALCULATIONS 
-CHECK  WHETHER  VALVE  FOUR  IS  CLOSED  OR  OPEN 
IF  < VALVS4.EG.0. )  GO  TO  340 
IF  (  DEL  !  A2  •  L,  I  .  CRACK 4  '  GO  TO  340 
-VALVE  FOUR  IS  OPEN  -  FLOU  RATE  CALCULATION 
-CHECK  FOR  CHOKED  FLOU  CONDITION 
IF  (VMACH2.GF-.!.)  GO  TO  330 
f F  ( DELTA2 . GT . CHANG4 )  GO  TO  510 
D°/-D1 
D6=D2 
GO  TO  320 
U5=D3 
D6=D4 

ORA  TE4=EXP ( D5T D6KAL0G (BELT A2 ) > 
0RATE4=GRATE4/60 .  *VAL  VS4*GC0RR2 
-CHOKED  VALVE  -  FLOU  RATE  CALCULATION 
RATEV4=GRATE4*DENS2 
GO  TO  400 

-VALVE  FOUR  IS  CLOSED 
RA!£V4=0.0 

-TOTAL  VALVE  FLOU  RATE 
RA  TE2=RATEV3+RATE V4 
FILTER  TWO  FLOU  RATE  CALCULATION 
-CHECK  FOR  CHOKED  FLOU 

IF  <  F  IL  IS2. EG . 0 . >  GO  TO  610 
IF  (VMACH3.GE. t . )  GO  TO  510 
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«RArF2=Fl+F2*RELTA3+F3*DELTA3**2+E  4*D£L!A3#*:J 

IF  (KEY3.EG.  1  )  QRATF2=0.00 

OF;  A 7 E2=QRATF2*AFIL  T2*FILT S2/60 . *QCORK2 

C - CHOKED  FILTER  -  FLOW  RATE  CALCULATION 

SiO  RATEF2=GRATE2*DENS2 

C=="*==ORIFIi:£  TWO  FLOW  CALCULAI'QN  -  ADD  IN  LFAK  CONT KI DU  1 1 ON  IF  ANY 

!- - CHECK  for  choked  flow 

610  IF  <aLEAKL'.E0.0.  )  GO  TO  700 
IF  VMACH3 .  GE  .  1 .  )  GO  10  620 

RATEL2=DC0EFF*ALEA|\2*S0RF  64 . 348*D£NS24 DEL T A3* 1 44  •  ) 

ORA fL2=RATEL2/DENS2 
GO  TO  700 

L - CHOKED  CONDITION 

620  RATEL2=QRATL2*DENS2 

l, - CALCULATE  THE  TOTAL  FLOW  RATE  AND  PRESSURE  DERIVATIVE 

700  RATF3=RAT£2+RATEL2+RA7£F2 
800  PREDER (2 ) =-MULT2#RAT£3 

c **#♦***####**##**#**♦#***♦♦*#**#♦♦##****##<****#♦♦****#*#*# t****##*Mtc 

C***M  OUTER  VOLUME  CALCULATIONS  **#*#C 

ctuntimmmHH  *********************  **t*M#**ts*Mt***Mt*t*t*«**c 

L - CALCULATE  GAS  PROPERTIES 

MASS 1=PR£SS1 /MULT 1 
HENS  1  -"MASS  1  /VOL  1 

C - CALCULATE  THE  ORATE  CORRECTION  FAC  I  OR 

0C0RR1=SQRT (PREXT( 1 ) /PRESS 1 ) 

OCORR 1 =S0RT ( QC0RR1 ) 

C - CALCULATE  THE  PRESSURE  DIFFERENCE 

KEY 1=0 

DELTAI=(PRESS1-PEXT)/144. 

IF  (DELTA1.LE.0. >  KEY 1=1 

IF  ( DEL  TA1 »LE  » 0. )  DELTA 1 =0 . OOOOOO 

C - CALCULATE  PRESSURE  RATIO  AND  THROAT  VELOCITY 

PRAT  1=F'£X  f /PRESS  1 
PEX!=PEXT/144. 

VET  1=1 .-PRATHurCXPl 

TF  (VEL1  .LE.O.  .*  VEL  1=0. OOOOOO 

VEL  1=SQF:T  ( VMUL ? *PRKSS1*VEL  1  /DENS1  ) 

IF  (VEL 1 .GT.VSOUND)  VEL1=VS0UND 
VMACHt-VEL 1 /VSOUND 
C=  =  --s-VALVE  ONE  FLOW  CALCULATIONS 

C - CHECK  WHETHER  VALVE  ONE  IS  OPEN  OK  CLOSED 

IE  VALVS1  .EQ.O.  )  GO  10  1040 
IF  (DEL TAi .LT. CRACK! )  GO  TO  1040 

. - VALVE  ONE  IS  OPEN  -  FLOW  RATE  CALCULATION 

C - CHECK  FOP  CHOKED  FLOW  CONDI  I  ION 

IF  (VMACH1 .GE. I . >  GO  TO  1030 
IE  (DFLIA1.G1.CHANGI !  GO  TO  1010 
AS  -A  1 
A6=A2 

GO  10  1020 
1010  A5=A3 

A6=A4 

1020  ORA TE 1 =FXP ( A5+Asir  ALOG ( DEL  !A!  )  > 

QPA I  £  1  =GK A  :  FT  1  760 .  * VAL  VS  1  *GCOPR  1 

f, - CHOKED  VALVE  -  FLOW  RATE  CALCULATION 

1 030  RA I EV 1 =QKA I  El *DENS 1 

GO  TO  1100 

L - VALVE  ONE  IS  CLOSED 

1040  RAT  EV1 =0 . 0 
C=====VALVE  TWO  El  OW  CALCOLATTONS 
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c - CHECK  WHETHER  VALVE  TWO  IS  CLOSED  OR  OPEN 

1100  IF  (VALVS2.EQ.0-)  GO  TO  1140 

IF  (DELTA1 .LT.CRACK2)  GO  TO  1140 

c - VALVE  TWO  IS  OPEN  -  FLOW  RATE  CALCULATION 

C - CHECK  FOR  CHOKED  FLOW  CONDITION 

IF  ( VMACH1 . GE . 1 . >  GO  TO  1130 
IF  ( DELTA 1 .GI .CHANG2)  GO  TO  1110 
B5=B1 
K6= 62 

GO  10  1120 
1110  Bb=B3 
B6~B4 

1 120  QRATE2-=EXP(B5+B6*AL0G (DELTA  1 )  > 

QRAT  E2=QRA TE2/60 . * VAL  VS2*UC0RR1 

C - CHOKED  VALVE  -  FLOW  RATE  CALCULATION 

1130  RATEV2=QRATE2*DENS1 
GO  10  1200 

r-— — VA1  Vt  TWO  IS  CLOSED 
1140  PA TEV2=0 . 0 

C===  =“F 1 1  TER  ONE  FLOW  RATE  CALCULATION 

C - (  HFCK  FOR  CHOKED  FLOW 

1200  >F  iFILTSl .EQ.O. )  GO  10  1410 
IF  (VKACH1.GE.1.)  GO  10  1310 

(1PATF1=E1  +E2*DELTA1+E3*DEL  TA  l  **2+E4*DELT»11*t3 

IF  (KEY1.E0.1)  ORATF 1=0 . 00 

(IRA  1F1=URAT  FlfAFILT  1*F  ILTS1/60.  #l«CORRl 

L - CHOKED  FILTER  -  FLOW  RATE  CALCULATION 

1310  KA TEF1 =QRATF 1 *DENS1 

C====-URIFIC£  ONE  FLOW  CALCULATION  -  ADD  IN  LEAK  CONTRIBUTION  IF  ANY 

C - CHECK  FOR  CHOKED  FLOW 

1410  IF  (ALEAK1 .EQ.O. >  GO  10  1500 
IF  (VMACH1.GE.1. >  GO  10  1420 

RA  TEL1  =  DC0EFF*ALEAK1»SQRT( 64 . 348*DENS1*DEL  TA1*144 . ) 

QRA IL1=RA TEL 1 /DENS 1 
GO  TO  1500 

C - CHOKED  CONDITION 

1420  RATEL1=«RATL1*DENS1 

C - CALCULATE  THE  TOTAL  FLOW  RATE  AND  PRESSURE  DERIVATIVE 

1 bOO  RATE 1 =RA  TE V 1 +RA TE V2+RATEL 1+RATEF 1 

PR£D£R( 1 >=-HULT 1 t ( RATE1-RATE2) 

RE  TURN 
END 
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Appendix  D 

Sample  Input  and  Output  Data  files 


Format  for  Input  Data  File  FTN30.DAT 


Line  Format 

1  20A4 

2  5A4 

3  6F10. 0 

4  20A4 

5  7  F 1 0.  0 


6  20A4 

7  7F10. 0 

8  20A4 

9  6F10. 0 


10  20A4 

11  F  10.0 


Title 
Gas  type 

Volume  #1,  volume  #2,  initial  pressure,  initial 
temperature,  gas  constant,  ratio  of  specific  heats 
Valve  type  #1 

Number  of  valves,  cracking  pressure,  knee  pressure, 
curve  coefficients  A  &  B  (below  knee  pressure), 
curve  coefficients  C  &  D  (above  knee  pressure) 

Valve  type  §2 

Same  as  line  5  except  for  valve  type  tt 2 
Filter  type 
Number  of  filters, 

filter  area/multiplier  (for  RA-2500  is  exit  area; 

for  CW-19  is  length  multiplier) 
curve  coefficients  A,  B,  C,  &  D 
Leak/orifice  title 
Effective  area  of  leak  or  orifice 
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12  4F10.0 

13  110 

14  8F10.  0 


Ending  time  of  calculations, 
timestep  =  0.  01  seconds, 
print  interval  multiplier  =  100. , 
accuracy  requirement  =  0.001 
Number  of  entries  to  follow  in  external  pressure 
table 

Time,  external  pressure 


Notes:  For  valves,  filters,  and  orifices  not  present  include  a 

title  but  leave  numerical  data  lino  blank. 

In  line  12,  the  timestep  and  print  interval  multiplier. 

determine  printout  time:  100  X  0.  01  =  1  second  printout. 

In  line  14,  4  pairs  per  line,  repeating  if  necessary. 

For  PRESSM.FOR,  repeat  lines  4  through  11  for  second 
volume  with  primary  volume  data  first  and  secondary 
volume  data  following  (see  IRBS  Payload  input  file). 

Units  for  Input  Data  File 

Gas  volume 
Initial  pressure 
Initial  temperature 
Gas  constant 

Valve  crack:ng  &  knee  pressure 
Filter  area 
Leak  area 
End  time 

External  pressure  t-ble: 
time 
pressure 
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cubic  feet  | 

pounds  per  square  inch  j 

degrees  F 

foot-pounds  per  pound-degree  R 

pounds  per  square  inch  j 

§ 

square  inches  I 

square  feet  | 

seconds  | 

seconds  | 

pounds  per  square  foot  1 


IP . . . . . . . * . . . . . . . . . . . . ■ . . . . . . . . . . . . . . . 


Sample  Input  Data  File  (FTN30.DAT) 
for  ZIP  Payload  -  PRESS4.F0R 


IP  PA*'  JOU  JNIERNAI  PRESSURE  HISTORY 
N  t !  ROOt -J 


<-,-00 

0.12 

35.0 

55.20  1 

.  400 

f.  1  hTi  E 

SEAL  P 

-249  0.10  PSI  MARKER 

CRACKING 

PRESSURE 

6 . 

0 . 0387 

0.10  10 

.8789  4 

.  7952  0 . 

9767 

0.4956 

ClRtl  I 

St  Al  P 

7-637  0, 

50  PSI  MARK EH  CRACKING 

PRESSURE 

0. 

0.325 

0.59  12 

. '900  17 

.3978  3. 

8647 

0.4786 

Mill  PURE  KA 

1.2  MXCP;OrtE  IER  PORE 

SIZE 

;>  B  0 

.11045  -0 

.007017  2.018104 

PtKt-tl  r  St  AL  - 

-  NO  LEAK 

100. 

o.ot 

1  00 . 

O.001 

M  v  0 

1  827.  / 

•  .•.) 

1826  1 

2.0 

1831 .0 

5.0 

1812.  i 

i  70P , 

5.0 

1  78  1 .  A 

6 . 0 

1 7.63.5 

7 , 0 

1739.4 

P  ... 

!  7 1 1 .  ! 

9.0 

1678.6 

10.0 

1642. 1 

■1.0 

1601 , 4 

J5S7.  I 

17.0 

1509.0 

14.0 

! 457.6 

15.0 

1 403 . 2 

)  6  .  •! 

l 746 . 5 

1  7.0 

i/87.9 

18.0 

1 22  7.  7 

19.0 

1166.5 

,Ki  ,  , . 

* 104.6 

21.0 

1042.1 

32 .0 

4  7« .  3 

2.5-0 

916,7 

'<1.0 

854 . 4 

25.0 

792.7 

26 , 0 

731 .9 

>7.0 

672  3 

OH .  0 

6  14.0 

29 . 0 

557-5 

10.0 

£02.8 

31  .0 

450.4 

v.o 

401.1 

33.0 

155.4 

34.0 

.313  0 

35.0 

274.0 

3'> .  0 

238.  i 

17.0 

305.9 

38.0 

176.7 

io.O 

150.5 

40.'. 

1  .17,  5 

41.0 

106.7 

42.0 

88.8 

43.0 

73.2 

,1  a  ,  i\ 

59.8 

45.0 

48,5 

46.0 

39,0 

47.0 

31,2 

48.0 

24.8 

49.0 

19.6 

50 . 0 

15.5 

51.0 

12. 1 

57-0 

<•'.4 

S3 . 0 

7.3 

54,0 

5.6 

55.0 

4.3 

56.0 

7.7 

5  7.0 

2.5 

58.0 

1 . 8 

59 . 0 

1.4 

AO .  0 

1.0 

61,0 

0.7 

62.0 

0. 5 

63.0 

0.3 

64.0 

0.2 

65.0 

0.2 

66.0 

0. 1 

67.0 

2,  i 

68.0 

0 .  Of. 

69.0 

0.0 

70.0 

u .  0 

360.0 

0.0 

Sample  Input  Data  File  (FTN30.DAT) 
for  1RBS  Payload— PRESSM. FOR 


IRPS  PAYLOAD  INTERNAL  PRESSURE  HISTORY 
AIR 

46.80  0.1  7 >7  0.50  70.0  53.35  1.400 

CIRCLE  SEAL  R7-637  tl.50  PSI  MARKED  CRACKING  PRESSURE 

3.  0.325  0.59  12.7900  17.3978  3.8647  0.4786 

NO  SECOND  VALUE  TYPE  PRESENT 

NO  FILTERS  PRESENT 

PERFECT  SEAL  —  NO  LEAKS 

CIRCLE  SEAL  P-249  0,10  PSI  MARKED  CRACKING  PRESSURE 

?,  0.0387  0.10  10.879°  4.7952  0.9767  0.4956 

NO  SECOND  VALUE  PRESENT 

NO  FIL  TF.RR  PRESENT  » 

DOOR  LEAK  DUE  TO  SEAM  SEAL 
0.000042 


IOC. 

89 

0.01 

100. 

0.001 

0.0 

182  7 , 7 

1.0 

1827.7 

7.0 

1816.7 

3.0 

1809.8 

4.0 

1798.8 

5.0 

1785.4 

6.0 

1768,2 

7.0 

1747.7 

9.0 

1723.2 

9.0 

1695.5 

10.0 

1663.9 

11.0 

1628.7 

12.0 

1589.6 

13.0 

1547,9 

14.0 

1502.1 

15.0 

1453.9 

16.0 

1 403-5 

17.0 

1348-8 

18.0 

1292.8 

19.0 

1-236.0 

20.0 

1178.8 

?  1 . 0 

1121.3 

22.0 

1063.6 

23.0 

1005.5 

24.0 

04'’ .  ^>5 

25 . 0 

889.07 

26.0 

830.89 

27.0 

773.17 

28.0 

715,94 

29.0 

659 . 67 

30.0 

604.40 

31.0 

550.43 

32 . 0 

498-10 

33 . 0 

447.64 

34.0 

400.17 

35.0 

355.76 

36.0 

3M.45 

37.0 

276 . 25 

33.0 

241.40 

39.0 

209.09 

40.0 

180.08 

41.0 

153.95 

42,0 

130.62 

43.0 

109.96 

44.0 

91 .890 

45.0 

76.198 

46.0 

62.720 

47.0 

51.207 

48.0 

41.468 

49,0 

33-300 

50.0 

26.506 

51.0 

20.907 

52,0 

16.354 

53.0 

12.708 

54,0 

9.8087 

55.0 

7.5328 

56.0 

5.7511 

57.0 

4.3639 

58.0 

3.2927 

59.0 

2.4401 

60.0 

1.6404 

61.0 

1.3554 

62.0 

0.9910 

63.0 

0.7203 

64,0 

0.5213 

65.0 

0.3756 

66.0 

0,2685 

67.0 

0.1897 

68.0 

0.1323 

69  0 

0.0°11 

70.0 

0.0617 

71.0 

0.0411 

72 . 0 

0.0269 

73 . 0 

0.0173 

74.0 

0.0111 

75.0 

0.0072 

76.0 

0 . 0046 

77.0 

0.0030 

78.0 

0.0020 

79.0 

0.0013 

80.0 

0.0009 

91  0 

0.0001 

32.0 

0.0000 

90.0 

0.0000 

100.0 

360-0 

0 . 0000 

0 . 0000 

150.0 

0 . 0000 

200.0 

0.0000 

350.0 

0.0000 
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Sample  Output  Data  File  (FTN31.DAT) 
for  ZIP  Payload  -  PRESS4.FOR 


ZIP  PAYLOAD  INTERNAL  PRESSURE  HISTORY 

gas  properties: 
type: 

MAIN  OOl UMfc 
INITIAL.  PRESSURE 
IFMPERATURE 
GAS  CONSTANT 


VALVE  uNfc  properties: 
type  : 

NUnBfP  fit  RELIEF  VALVES 
CRACK ! NO  PRESSURE 
CURVE  CHANGE  POINT 
CuEFF  tC LENT  1 
COEFFICIENT  2 
COEFFICIENT  3 
COfcFE  If. IEN 1  4 


NITROGEN 

6 .80  Cl  I  F! 

0 . J  2  PSI 
35.0  DEGREES  F 
55.20  f I-L R/L b-DEG  R 


CIRCLE  SEAT  R-J49  010  PSI  NARKED 

6*  I'RACMNSi  PRESSURE 

O.o4  FSL 
0.10  r’S  ( 

10  8?9 
4.  ?95 

4.D7? 

0 . 4^6 


VALVE  LUO  properties: 
type: 


CIRCLE  SEAL  R7-657  v'  50  Pfc>1  HANKED 


FILTER  properties: 

type:  MTCl-tPORE  RA  1.2  MICROMETER  PORE  SIZE 
NUMBER  OF  FILTERS  =  CRACKING  PRESSURE 
EXIT  AREA  =  0.11045  SO  IN 


COEFFICIENT 

1 

=r 

-0.G0‘7 

COEFFICIENT 

0 

= 

2.018 

COEFFICIENT 

3 

•r 

0.000 

COEFFICIENT 

4 

= 

0.000 

POOR  LEAK  PROPERTIES: 

type  of  seal: 

PERFECT  SEAL.  —  NO  L  EAK 

choking  properties: 

RATIO  OF  SPECIFIC 

HEATS 

=  1.400 

CRITICAL  PRESSURE 

RATIO 

-  0.5283 

SF’EED  OF  SOUND 

1012.4  FI’S 

TTMF 

EXTERNAL 

PRESSURE 

INTERNAL 

RRESSURE 

DEI TA  P 

INTERNAL 
GAS  MASS 

TUT  A!  MASS 
Fl  OM  KALE 

SECS 

PSI 

PSI 

PSI 

LBM 

l  BM/SEC 

0.0 

12.6V 

12.81 

0.120 

0,45946 

0 . 006  1326 

1.0 

12.68 

12-75 

0.066 

0.45/12 

a . 0008095 

2.0 

12.65 

12.  72 

0.072 

0 . 45606 

0.0012052 

3.0 

12.59 

12.67 

0 . 08 1 

0.45423 

0  -  0021200 

4  0 

J2.50 

12.5? 

0.088 

0.45136 

0 . 0030970 

5.0 

T2.3V 

12.48 

0 . 093 

0. 44751 

0.0040261 

6.0 

12.25 

12.34 

0 . 098 

,  44267 

0- 0049980 

7.0 

12.08 

12. 1? 

0.108 

>' ,  43702 

0  0057598 

69.0 

0.00 

0.33 

0.333 

0.01196 

0.0008620 

69  .  0 

0 , 00 

0.31 

0.310 

0.01113 

0.0008020 

70  0 

0.00 

0.39 

0.289 

0.01036 

0.000  7-163 

71  -0 

0 . 00 

0.27 

0.269 

0.00964 

0.0006944 

/:> .  o 

0.00 

0.25 

0 . 250 

0.00897 

0.0006461 

7,5.0 

0.00 

0.23 

0.233 

0.00834 

0.0006012 

74.0 

0.00 

0.22 

0.216 

0.00776 

0.0005594 

/S .  /) 

0 . 00 

0.20 

0.201 

0.00722 

0.0005205 

//,  ,  y 

0,00 

0.19 

0.187 

0.00672 

0.0004843 

7>,0 

0.00 

0.1? 

0.174 

0.00625 

0.000450? 

78,0 

0 . 00 

0.16 

0.162 

0.0O582 

0.0004193 

79 . 0 

0.00 

0.15 

0.151 

0.00542 

0.0003902 

no ,  0 

0.00 

0,14 

0.141 

0.00504 

0.0003631 

81.0 

0 , 00 

0.13 

0.1  51 

0.00.469 

0.0003378 

8  .’.0 

0 . 00 

0.12 

0.122 

0.00436 

0.0003143 

87 ,  v) 

0 . 00 

0.11 

0,113 

0.00406 

0.0002925 

8-1.0 

0.00 

0.11 

0.105 

0.00378 

0.0002721 

85,0 

0.00 

o.  io 

0.098 

0.00351 

0.0002532 

86.0 

0.00 

0.09 

0.091 

0.00327 

0 . 0002356 

8 7  0 

o.oo 

0.08 

0.085 

0.00304 

0.0002192 

88  0 

0 . 00 

0,08 

0.079 

0.00283 

0.0002040 

89 . 0 

o .  oo 

0.0? 

0.073 

0.00263 

0.0001899 

90-0 

o .  00 

0.07 

0 . 068 

0.00245 

0.0001766 

91.0 

0.00 

0.06 

0.064 

0.00228 

0.0001643 

9'.’ .  0 

0-00 

0.06 

0.059 

0.00212 

0.0001529 

93 .  0 

0.00 

0.06 

0.055 

0.00197 

0.0C01423 

9-1.0 

0.00 

0  •  05 

0.051 

0.00184 

0,0001324 

95. 0 

0.00 

0.05 

0.048 

0.00171 

0.0001232 

96.0 

0.00 

0.04 

0  "44 

0.00159 

0.0001146 

97.0 

0.00 

0.04 

0.00148 

0.0001067 

98.0 

0.00 

0.04 

0.  „9 

0.00139 

0.0000035 

99.0 

0.00 

0.04 

0.039 

0.00138 

0.0000035 

100.0 

0.00 

0.04 

0.038 

0.00138 

0.0000035 
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A 

a 

C 

c 

P 

M 

m 

m 

n 

P 

pc 

Q 

R 

T 

t 

V 

v 

y 

p 


Nomenclature 

Area 

Speed  of  sound 
Constant 

Pressure  coefficient 
Mach  number 
Mass 

Mass  flow  rate 
Correction  factor  exponent 
Pressure 

Cracking  pressure  differential 

Volume  flow  rate 

Gas  constant 

Temperature 

Time 

Volume 

Velocity 

Ratio  of  specific  heats 
Density 


Subscripts 

atm  Atmospheric  conditions 

crit  Critical  (at  M=  1) 

e  Exit,  external 

f  Filters 

i  Inlet,  internal 

1  Leaks 

o  Orifices 

T  Total 

v  Valves 

1  Primary'  volume  internal 

2  Secondary  volume  internal 
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